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■ Properties of dense nucleon matter and the structure of neutron stars are 
studied using variational chain summation methods and the new Argonne v is 
two-nucleon interaction, which provides an excellent fit to all of the nucleon- 
nucleon scattering data in the Nijmegen data base. The neutron star gravita- 

■ tional mass limit obtained with this interaction is 1.67M & . Boost corrections 
CN ■ to the two-nucleon interaction, which give the leading relativistic effect of or- 
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and 2.20 M & . Hamiltonians including a three- nucleon interaction predict a 
transition in neutron star matter to a phase with neutral pion condensation 
at a baryon number density of ~ 0.2 fm -3 . Neutron stars predicted by these 
Hamiltonians have a layer with a thickness on the order of tens of meters, over 
which the density changes rapidly from that of the normal to the condensed 
phase. The material in this thin layer is a mixture of the two phases. We 
also investigate the possibility of dense nucleon matter having an admixture 
^ ■ of quark matter, described using the bag model equation of state. Neutron 

stars of IAMq do not appear to have quark matter admixtures in their cores. 
However, the heaviest stars are predicted to have cores consisting of a quark 
and nucleon matter mixture. These admixtures reduce the maximum mass 
of neutron stars from 2.20 to 2.02 (1.91) M for bag constant B = 200 (122) 
MeV/fm 3 . Stars with pure quark matter in their cores are found to be unsta- 
ble. We also consider the possibility that matter is maximally incompressible 
above an assumed density, and show that realistic models of nuclear forces 
limit the maximum mass of neutron stars to be below 2.5M . The effects 
of the phase transitions on the composition of neutron star matter and its 
adiabatic index V are discussed. 

PACS numbers:21.65.+f, 26.60.+C, 97.60.Jd 
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I. INTRODUCTION 



The significant influence of nuclear forces on neutron star structure is by now firmly 
established by a large body of theoretical and observational evidence [|l|. In the absence 
of these forces, the maximum possible mass of neutron stars composed of non-interacting 
neutrons is ~ 0.7 solar masses (M ) p[. Since most observed neutron star masses are above 
1.3M [Q, they must be supported against gravitational collapse by pressure originating 
from nuclear forces. In the present work, we study neutron star structure using one of 
the most realistic models of nuclear forces currently available. A brief outline of previous 
calculations leading to this work is presented below. 

Shortly after the discovery of pulsars, calculations of the equation of state (EOS) of 
neutron star matter with realistic models of the two nucleon interaction (NNI), obtained by 
fitting the nucleon-nucleon (NN) scattering data then available, were carried out using the 
lowest order constrained variational method jH|. The results demonstrated that nuclear 
forces increase the mass limit of stable neutron stars beyond 1.4M . 

By the late 1970s it had become clear that the NNI alone could not account for the 
properties of nuclear matter or few-body nuclei. Variational P and Brueckner calculations 
@, including higher order cluster contributions, established that nuclear matter with realistic 
NNI saturates at too high a density. In addition, these interactions were known to underbind 
3 H. Ignoring the latter problem, plausible density dependent terms were added to the Urbana 
v 14 (U14) model of NNI |S||J to reproduce the observed equilibrium properties of nuclear 
matter. The resulting density dependent (U14-DDI) model of nuclear forces predicted stable 
neutron stars having masses up to 1.8M HlQ| , [iT| . 

Since nucleons are made up of quarks and have internal degrees of freedom, we can 
expect interactions among three (and perhaps four or more) nucleons, in addition to the 
NNI. The Urbana three nucleon interaction (TNI) models contain only two terms, with 
strengths fixed by the saturation density of nuclear matter and the binding energy of 3 H. 
Wiringa, Fiks and Fabrocini (WFF) used the U14 and the subsequent Argonne 
(A14) 0] models of NNI, together with the Urbana VII (UVII) model of TNI, to study 
neutron star structure and obtained mass limits of 2.19M and 2.13M with the U14+UVII 
and A14+UVII, respectively. They also found that pure neutron matter (PNM) undergoes 
a transition to a phase having spin-isospin order, attributed to neutral pion condensation, at 
a density of ~ 0.2 fm" 3 with the A14+UVII, but not with the U14+UVII. Neither of these 
models results in a phase transition in symmetric nuclear matter (SNM), which is composed 
of equal numbers of neutrons and protons. 

In the early 1990s the Nijmegen group |14| examined carefully all NN scattering data at 
energies below 350 MeV published between 1955 and 1992. They extracted 1787 proton- 
proton (pp) and 2514 proton- neutron (np) "reliable" data, and demonstrated that these data 
determine all NN phase shifts and mixing parameters quite accurately. The NNI models 
which fit this Nijmegen data base with a x 2 l^data ~ 1 are called "modern". These include 
the Nijmegen models [T5| : Nijmegen I, II and Reid-93, the Argonne Vi 8 [JTH], denoted here 
by A18, and the CD-Bonn fIT|| . In order to fit the pp and np data simultaneously and 
accurately, these models include a detailed description of the electromagnetic interactions 
and terms that break the isospin symmetry of nuclear forces. All include the long range 
one-pion exchange potential, but follow different treatments of the intermediate and short 
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range parts of the NNI. The differences among the predictions of these models for the 
properties of many-body systems are much smaller than those among the predictions of 
older models, presumably because all modern potentials accurately fit the same scattering 
data. For example, the 3 H binding energies predicted by the modern Nijmegen models and 
A18 are between -7.62 to -7.72 MeV @, while that of CD-Bonn is -8.00 MeV flTj. The 



difference between these results and the experimental value of -8.48 MeV is used to fix one 
of the parameters of Urbana TNI models. 

Detailed studies of the energies of dense nucleon matter were carried out recently by 
Engvik et al. \V5\ using all the modern models of NNI and the lowest order Brueckner 



(LOB) method. According to these studies, the results with the modern potentials are all 
quite similar up to densities relevant to neutron stars. For example, the energies predicted 
with the LOB method for the energy of neutron matter at 5po, where po = 0.16 fm -3 is the 
equilibrium (saturation) density of nuclear matter, range from 80 to 93 MeV per nucleon. 
The spread of 13 MeV in these energies is small compared to the possible errors in the LOB 
method and the expected contributions of TNI at this density. This model independence 
results from the fact that the mean interparticle distance at p ~ 5po is greater than 1 fm, 
and the predicted matter energy is therefore not sensitive to the details of the interaction 
at < 1 fm. 

In this paper we study the structure of neutron stars with the A18 model using variational 
chain summation (VCS) methods, which hopefully include all leading many-body correlation 
effects. The Urbana model IX (UIX) [p(J is used to estimate the effect of TNI. Previous 



studies of nucleon matter with A18 and UIX interactions have indicated the possibility of a 
transition to a neutral pion condensed phase for both PNM and SNM | 2~I|| . The effects of 
such a transition on the structure of neutron stars are studied here in detail. The effect of 
relativistic boost corrections [2~2~| to the A18 interaction is also examined. At high densities, 
we consider the possibility of matter becoming maximally incompressible, as well as that of 
a transition to mixed phases of quark and nucleon matter 



The relativistic mean field (RMF) approximation |?5] has been used in many studies of 
high density matter and neutron stars. There exists a vast amount of literature on this topic, 



some of which has been reviewed by Glendenning [26|. While the RMF approximation is very 



elegant and pedagogically useful, it is not valid in the context of what is known about nuclear 
forces, which is the theme of this work. For example, using the meson parameters of the 
CD-Bonn S-wave potentials in the RMF approximation leads to unbound SNM. (In order to 
accurately fit the NN scattering data, the phenomenological scalar meson parameters in the 
CD-Bonn model are allowed to depend on partial wave quantum numbers.) At a density of 
Po, for SNM, the RMF approximation yields an energy per nucleon of ~ +20 MeV, while the 
LOB method gives ~ — 18 MeV. In the mean field Hartree approximation, which is implicit 
in the RMF calculations, the A18 NNI gives energies per nucleon of ~ +30 (+37) MeV at po 
and +155 (+204) MeV at 5po fc> r SNM (PNM), while the variational calculations presented 
here give ~ -18 (+12) and +25 (+88) MeV, respectively. 

The main problem is that the mean field approximation for meson fields is only valid 
for pf <C 1, where p is the inverse Compton wave length of the meson and f is the mean 
interparticle spacing. Over the 1-5 po density range, f estimated using a body centered cubic 
lattice ranges from 2 to 1.2 fm. Thus pf is in the range 1.4 to 0.8 for the pion and 7.8 to 
4.7 for vector mesons. The mean field approximation is not applicable, since these values 
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are obviously far from being much smaller than one. The RMF approximation can be based 
on effective values of the coupling constants that take into account the correlation effects. 
However, these coupling constants then have a density dependence, and a microscopic theory 
is needed to calculate them. 

This paper is organized as follows. Section II contains a summary of the non-relativistic 
calculations with A18 and A18+UIX models of nuclear forces, while section III describes 
the calculations including the relativistic boost interaction, denoted by 5v, without and with 
the TNI model UIX*. The beta equilibrium of neutron star matter is discussed in section IV 
and results for neutron star structure are presented in section V, where we also discuss the 
effects of the possible transition to mixed nucleon and quark matter phases. The adiabatic 
index and sound velocities in neutron star matter are given in section VI, and conclusions 
are presented in section VII. 



II. NONRELATIVISTIC CALCULATIONS 

Nonrelativistic calculations of SNM and PNM with the A18 and UIX interactions were 
carried out using Variational Chain Summation (VCS) techniques described in detail in 
||2 1| | . Energies are calculated by evaluating the expectation value of the Urbana-Argonne 
Hamiltonian with a variational wavefunction composed of a product of pair correlation op- 
erators acting on a Fermi gas wavefunction. The pair correlation operators are written as 
a sum of eight radial correlation functions, each multiplied by one of the two-body opera- 
tors: (1, (Ti ■ a-j, Sij, (L ■ S)ij) ® (1, Ti ■ Tj). The wavefunction depends on three variational 
parameters: the range of the tensor correlations, d t , the range of all other correlations, d c 
and a quenching parameter a, meant to simulate medium effects. In this section we discuss 
results obtained for four cases, namely SNM and PNM with and without the three-nucleon 
interaction, in order to indicate their sensitivity to various terms in the nuclear force, and 



to extend them to higher densities, beyond the range covered in [^T[ . 

The optimum values of the parameters d t and d c in matter without and with the three- 
nucleon interaction are shown in Fig. [|. Some of the noise in the variation of these parameters 
with matter density is due to the insensitivity of the energy to their values at the variational 
minimum. The large increase in dt of PNM without Vijk at p ~ 0.5 fm~ 3 is due to a 
transition to a phase with 7r°-condensation as discussed in [^TJ. With V^k there are sudden 
changes in the d c and d t of PNM at p ~ 0.2 fm~ 3 , and in the d t of SNM at p ~ 0.32 fm" 3 . 
These are associated with the same phase transition. We note that the present variational 
wavefunction is not fully adequate to describe the long range order in the 7r°-condensed 
phases. However, since the change in the energy due to the 7r°-condensation appears to be 
small, we expect our estimates of the energies of matter with 7r°-condensation to be useful. 

Our nonrelativistic Hamiltonian, Hnr, is comprised of the nonrelativisic kinetic energies 
and the two-body A18 and three-body UIX interactions. The NNI includes a static, long- 
range one-pion exchange part, with short-range cutoff, and phenomenological intermediate 
and short-range parts, which depend on the six static two-body operators (1, <Xj • (Tj, Sij) ® 
(l,Tj ■ Tj) and eight momentum-dependent (MD) two-body operators: ((L • S)^, Lf,-, <Ti ■ 
(TjLfj, (L • S)h) ® (1, Ti ■ Tj). The A18 includes an additional isovector operator, and three 
isotensor operators, which distinguish between pp, np and nn interactions. These isovector 
and isotensor terms are small, and give zero contribution to the energy of SNM to first 
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order. They are therefore neglected in SNM calculations. In the case of PNM all the isospin 
operators can be eliminated and the full A18 with isovector and isotensor terms becomes 
the sum of a static part with operators (l,(Tj • (Tj,Sij), and a MD part with operators 
(L • S)ij, Lip o~i ■ CjLfj and (L • S)?-. The UIX model of contains two static terms; the 
two-pion exchange Fujita-Miyazawa interaction, V^., and a phenomenological, intermediate 
range repulsion V-%. The strength of the V^. interaction was determined by reproducing 
the binding energy of the triton via Greens Function Monte Carlo (GFMC) calculations pU| . 
while that of V^ k was adjusted to reproduce the saturation density of SNM. 

Expectation values of the various interactions are calculated in the VCS framework 
by summing terms in their cluster expansions. The one-body, two-body and many-body 
contributions to the kinetic and NNI energies are listed in Tables FHlVl- The one-body 
cluster contribution includes only the Fermi gas kinetic energy, Tp. The remainder of the 
kinetic energy is separated into the contribution from the two-body cluster {T) 2 b and that 
coming from the many-body clusters, (T) MB . The kinetic energy can be calculated using 
different expressions related by integration by parts. If all MB contributions are calculated, 
these expressions yield the same result. However, they yield different results when only 
selected parts of the MB clusters are summed by VCS techniques. We have calculated the 
many-body kinetic energy using expressions due to Pandharipande and Bethe (PB), and to 
Jackson and Feenberg ( JF) . The averages of the PB and JF results appear under the 
column (T) MB , and the differences T PB — T JF are listed under AT. In studies of atomic 
helium liquids the exact energies, calculated via Monte Carlo (MC) methods, lie between 



the PB and JF values evaluated using VCS methods p8 |. Although exact MC calculations 



are as yet not practical for nucleon matter, we believe the average of the two expressions 
to be more accurate than either one, and that the difference provides a measure of the 
uncertainty in the many-body calculation. The (T)mb is quite small in nucleon matter, 
due to cancellations between various many-body terms, and therefore the difference AT is 
a better indication of this uncertainty. 

The two-body cluster contribution to the static and MD parts of the NNI energy are listed 
under (v s ) 2 b, and (v md ) 2 b- The (v s ) 2 b is negative and large enough to bind SNM, though 
not PNM. The (v md ) 2 b increases rapidly with density, and is proportional to pT F at p > 2p . 
The many-body contributions to the static and MD parts of the NNI energy are listed 
separately as (v s )mb and (v md )mb- Previously, great efforts were made to improve upon 
the accuracy of the calculation of (v s )mb HHf; however, at large densities the (v md )mb 
grows rapidly, and becomes larger in magnitude than (v s )mb- The MD contribution is more 
difficult to calculate because of the gradients in v^ D , which may operate on the correlations 
of nucleons i and j with other nucleons. All the leading terms are calculated as discussed 
in [^TJ, and the corresponding errors should therefore be much smaller than the reported 
values. At higher densities the (v md )mb becomes proportional to p(T) 2 B, as expected. 
An additional perturbative correction to the two-body energy is listed as 5E 2 b- This small 
correction is due to an improvement in the variational wavefunction, which occurs when 
correlation functions are calculated separately in each I, S, J channel pif . 

The expectation values of v^, v^, V?l and V^ k for SNM and PNM are listed in Table [V]. 
Here, v$ is the phenomenological part, — vfj, of the NNI. Since the UIX Vijk is purely 
static, the error in the calculation of its expectation value is likely to be small. In SNM, 
with only the A18 interaction, the v n gives more than half of the total NNI energy at all 
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densities considered. The corresponding calculation with A18+UIX interactions, shows a 
significant increase in the magnitudes of the negative contribution of VyjJ between p = 0.32 
and 0.4 fm -3 , associated with pion condensation. In PNM these pion exchange interactions 
make relatively small contributions at densities below the phase transition, occurring at 
p ~ 0.2 (0.5) fm~ 3 with (without) UIX. However, at the densities above the transition they 
make large negative contributions comparable to those in SNM. At the highest densities the 
contributions of V?fi. and V^ k become very large, and the validity of this purely nonrelativistic 
approach becomes questionable. As discussed in the following section, approximately 40 % 
of the contribution of V^ k is due to relativistic boost corrections to the NNI , and a more 
plausible theory is therefore obtained by removing this boost contribution from the UIX 
interaction. 

The total energies, calculated in the manner previously described for SNM and PNM, 
appear in Tables |VI| - |VII| , and in Figs. ||] and The pronounced kink in the E(p) of SNM 
with Vijk, at p = 0.32 fm -3 , is due to the phase transition; the corresponding feature in 
PNM is a somewhat more subtle change in the slope of the curve at p = 0.2 fm -3 . 

It is evident from the E(p) figures for SNM that without the Vijk, the present calculation 
cannot explain the empirical saturation density po, of nuclear matter. As previously noted, 
the strength of VM. is adjusted to obtain the correct equilibrium density in calculations 
with Vijk- However, the present calculations with V^k underbind SNM at saturation density, 
giving E(pq) ~ —12 MeV per nucleon instead of the empirical value of ~ —16 MeV. This 
discrepancy is presumably due to use of imperfect variational wavefunctions, which do not 
include, for example, three- and higher-body correlations. It is known from comparison 
of the results of variational Monte Carlo (VMC) and exact GFMC calculations pO| , |3T 



that variational wave functions of the present form underbind the light p-shell nuclei. The 
variational energy of 8 Be for example, is above the exact GFMC result by ~ 12 %, even 
after incorporating into the wave function some of the three-body correlations, which we 
have neglected here. 



III. RELATIVISTC BOOST CORRECTION TO THE NN INTERACTION 

In all analyses, the NN scattering data is reduced to the center of mass frame and 
fitted using phase shifts calculated from the NNI, v^, in that frame. The obtained by 
this procedure describes the NN interaction in that frame, in which the total momentum 
Pij — Pi +Pj3 is zero. In general, the interaction between particles depends upon their total 
momentum, and can be written as 

v{Pi J ) = v ij +5v{P ij ), (3.1) 

where i>y is the interaction for P^- = 0, and 5u(Py) is the boost interaction which is 
zero when Py = 0. 

It is useful to consider a familiar example. The Coulomb-Breit electromagnetic interac- 
tion |32| between two particles of mass m and charge Q, ignoring spin dependent terms for 
brevity, is given by 

^.p^SO-^-W)' (3 - 2) 
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up to terms quadratic in the velocities of the interacting particles. In our notation it is 
expressed as 



v (Pi, Pi) = v ij + $ v (Pij), (3-3) 

with 

^ St 1 < 3 - 4 ' 

5v(P tJ ) = - 5L -2L + 1 ; J - f , 3.5 

where = (p$ — Pj)/2 is the relative momentum. 

In all realistic models of v^, such as the A18, the dependence on is included in the 
momentum-dependent part of the interaction, v^f D . However, we have neglected the 5v(Pij) 
in the calculations presented in the previous section. Even though contributions of the boost 
interaction to the binding energy of SNM and 3 H were estimated by Coester and coworkers 
years ago |33|j34]] , these contributions have been neglected in most subsequent studies of 



dense matter. 

Following the work of Krajcik and Foldy p5| , Friar []3B| obtained the following equation 



relating the boost interaction of order P to the interaction in the center of mass frame: 

6v(P) = -pjo + ^[PtP-V,.] + [fa _ a ) x P ■ V, v}. (3.6) 
8m z 8m z 8m z 

The general validity of this equation in relativistic mechanics and field theory was recently 
discussed |p2|] . Incorporating the boost into the interaction yields a nonrelativistic Hamil- 
tonian of the form: 

h*nr = E + £K- + H?ij)) + E^3* + ---» ( 3 - 7 ) 

where the ellipsis denotes the three-body boost, and four and higher body interactions. 
This H^r contains all terms quadratic in the particle velocities, and is therefore suitable for 
complete studies in the nonrelativistic limit. 

Studies of light nuclei using the VMC method |37]|J^] find that the contribution of the 
two-body boost interaction to the energy is repulsive, with a magnitude which is 37% of the 
Vjj k contribution. The boost interaction thus accounts for a significant part of the V$ k in 
Hamiltonians which fit nuclear energies neglecting 5v. 

In the present calculations we keep only the terms of the boost interaction associated 
with the static part of Vij, and neglect the last term in Eq. ( |3.6|) . That term is responsible 
for Thomas precession and quantum contributions that are negligibly small here J39|. Our 
Sv is given by 

5v(P) = -^v s + -Lp • r P ■ W. (3.8) 
8m z 8m z 

The two terms are due to the relativistic energy expression and Lorentz contraction, and 
are denoted 5v RE and 5v LC , respectively. The three-nucleon interaction used in the H^ R 



7 



Eq. (|3.7|) is denoted by V£ k . Its parameters are obtained by fitting the binding energies of 
3 H and 4 He, and the equilibrium density of SNM, including 5v. The strength of VM is 0.63 
times that of V^ k in UIX, while that of V?j k is unchanged. The resulting model of Vij k is 
called I IX' . 

The approximate Hamiltonian Hnr, containing A18 and UIX interactions without Sv, 
and the more correct Hamiltonian H^ R , containing A18, UIX* and Sv interactions, yield 
very similar results for light nuclei up to 8 Be [3~8~| and for SNM up to equilibrium density. 
However, the two models differ at higher densities, since the contributions of Sv and V^ k 
have different density dependences. 

One may also consider relativistic nuclear Hamiltonians of the type 



Hr = J2 sjvl + ™ 2 + E(%' + 8v(?ii)) + E Vijk + ■ ■ ■ , (3-9) 

which require re-fitting the two-nucleon scattering data to determine the two-body interac- 
tion, Vij, using relativistic kinetic energies [|T7|| . In light nuclei, the Sv contribution accounts 
for most of the difference between the energies obtained with H R and H NR , since the dif- 
ference between the contributions of the nonrelativistic and relativistic kinetic energies is 
largely cancelled by the difference in interaction energy contributions from and v^. The 
results obtained with H^ R are very close to those from H R , indicating that the former 
represents a significant improvement over H^r 



A. Calculation of Boost Interaction Energy 

The relativistic boost contributions are calculated by evaluating terms in the cluster 
expansion of {Sv^ E ) and (Svfj ). In addition to the dominant two-body cluster, we have 
calculated dressed three-body separable diagrams and central chain diagrams. 

In the case of the two-body cluster, the gradients in the center of mass momentum 
operator, Pjj = — z(Vj + Vj), can act only on the Fermi gas part of the wavefunction, since 
the correlations depend only on the relative coordinate. Thus, the two-body cluster 
contribution to (Sv^ E ) is 

= -^pE/ d% (/vr% c {oioTpi) (3-10) 

pmp' 

°™ {ex) = £ / rf % 0' 2 - M\ (f P v m f% C {omPTpti) ■ (3-11) 

pmp'n 

The quantities C(- • •) in the integrand of these expressions represent the spin-isospin- 
independent part (or "C-part") of the operator product enclosed by the parentheses. Only 
the C-part of operator products appear in the cluster integrals, since the energy expecta- 
tion value requires a sum over all possible a z and t z 0, which average to zero in isotropic 
matter. In the case of SNM the indices p,m,p' run over the first six operators, as we boost 
only the static interactions and consider only static correlations in this calculation. The 
index n=l,4 comes from the exchange operator. In the case of PNM the Tj • Tj operators 
are eliminated from the Hamiltonian, thus the indices p, m, p' and n can represent only unit 
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or spin-dependent operators. The exchange operator also contributes a factor of 1/s, where 
s is the degeneracy of the system (4 for SNM, 2 for PNM.) 

The C^(dir) differs from the corresponding expression for the two-body direct contri- 
bution to (vfj) by a factor — (l/8m 2 )(6£;f,/5), which is the expectation value of —P^/ (8m 2 ) 
in the Fermi gas. The exchange part of the cluster, Cj]f (ex), has the same form as the 
corresponding expression for (v^)2B(ex), with I 2 replaced by — (2)/(8m 2 )(/' 2 — /V 2 /), where 
/ = l(k F r) is the Slater function. This expression results from from gradients in acting 
on the plane waves 



2p \(l'(k F rij)) 2 - l(k Frij )V 2 l(k Frij )] . (3.12) 



An , 



Here, ^4 is the number of nucleons, Q is the normalization volume, and since we are in the 
thermodynamic limit, A/Q = p. 

The two-body cluster contribution to (Sv^°) is 

Ci£{dir) = *Z Y. I d% (rr—f) C (OJOffOj) (3.13) 

pmp' \ / ij 

C£ = £ / d% (I' 2 - ir) (pr^A C {0^0^ , (3.14) 

pmp'n \ / ij 



which is simply the cluster contribution of a nonrelativistic potential r dv/dr, with the 
direct term multiplied by (2/5)(k F /8m 2 ), and with the I 2 in the exchange term replaced by 
(2/8m 2 )(/ /2 — 11"). As in the 8v RE case, the extra factors in (Sv LC ) result from the gradients 
in Sv LC acting on the Fermi gas part of the wave function. 

The three-body separable diagrams represent the most significant many-body cluster 
contributions to the boost energy. The direct term for a boost 5v x , where x = RE or LC, 
has the form: 



£ E / (l '' r ' (l '' r 'i '/ :i /'V " k "' i ^ '' K H : 

■■I- / . _ / ** 



ijk pmp'qq' 



The plane waves are written here in terms of the relative momentum, = (k, — k,)/2, 
and the center of mass momentum, ICy = k« + k,, of the interacting pair. The interacting 
exchange (ij — ex) and the passive exchange (ik — ex) expressions are obtained from the 
above by inserting the appropriate exchange operators ((1/s) Z)n=i,4^- or (V s ) S n =i,4^«;) 
to the far left of each operator product, and replacing the first plane wave product by 

( -k, r, k r • K H ■ Qr ( '( k k 3 i- ,-K„ ir| 



Following the notation used in the calculation of the MD interaction energy |f40| , |2T 
separable diagrams are classified as K-diagrams and F-diagrams. The former have gradients 
in 5v x acting on the Fermi gas part of the wavefunction, and the latter have them acting 
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on the correlation operators F^. As with the two-body cluster contribution Cfg, the K- 
diagrams depend linearly on the Fermi kinetic energy. While Cf B scales as p Tp, like the 
(v MD )2B, the K-diagram contributions scale roughly as p 2 Tp. The K-diagrams generally 
make only small contributions to (Sv x ), the major separable contributions coming from F- 
diagrams. The relatively large contribution of the F-diagrams, versus the K-diagrams, can 
be understood in the following way. The correlated particle k in the separable diagram 
modifies the center of mass momentum of the interacting pair ij via Fik, thus enhancing 
the boost correction. As the form of the F-diagram integrals suggest, we find that their 
contributions exhibit the same scaling behavior as (v md )mb, namely as p(T) 2 B- 

K-diagram contributions to (Sv RE ) have been evaluated for the direct three-body sepa- 
rable diagram and the interacting exchange diagram. These contributions factorize into an 
integral over rij, which is simply the corresponding two-body boost diagram, and an integral 
over r ik . The latter integral is a so-called single- loop vertex correction, which is included in 
the more general vertex correction, M dd — 1, to (fl-), defined in |J. The direct K-diagram 
contribution is 



W S RE (K, dir) 



kp 3 



3m 2 5 



P 12 J d 3 r lJ ( K f p v m f p ')^K pmp ' A p ' 

A« P - J d 3 r tk {flf (D pq + D mq + D plq ) , (3.16) 



pmp'q ' 



where the K pmp , A p and D pq 
the non-commutativity of operators anu. 
contribution is given by 

W RE (K,ij -ex) = 
8m 2 s 



matrices, defined in ||, give the C-parts taking into account 
" and Oli.. The corresponding interacting exchange 



J . ^'2 _ /y 2 ^ ^ JP v m jp'^j \ ^j^npm' j^mp'm' _|_ j^prnm' j^-p'nm'^ j^m! 

P J d 3 r lk (ff k f A* (D qm , + l -D qm + ^D qn }j . (3.17) 



The F-diagram contributions to (5v RE ) are evaluated for the direct, interacting exchange 
and passive exchange terms. These contributions also factor into separate integrals, where 
the Tij integral has the form of the two-body contribution to (v s ), and the r ik integral is 
a new type of vertex correction involving gradients of the fik- The direct and r/'-exchange 
integrals are 



W s RE (F,dir)= J2 P -^d 3 r i3 {j p v m f p ').K pmp 'A p ' 



pmp'q 



Am 2 



J V T / ik 



- (D pq + D mq + D plq ) + 1 



A q 



(3.18) 



W RE (F, ij-ex)= ~— I (l 2 f P v m f p ') . . - [K npm 'K mp ' m ' + K pmm ' ' K p ' nm ') A m ' 



Am? 



p'inp' qm'n 

6 



/ d 3 r lk fl (Vf - (5 q5 + 5 q6 : 



tk 



2 D qm i -\- F) qm -\- F) qn -\- 1 



A q (3.19) 
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The i/c-exchange contribution is evaluated only to leading order, namely for the cases where 
at least one f ik is a central link. In this approximation, this contribution takes the form: 

W* E (F, ik-ex)= Yl -\\ [f v v m f) K pmp 'A p ' 

pmp'n 



Am 



1 + ~D np 



In this equation, the index n runs from 2 to 4 only. 



(3.20) 



The K-diagram separable three-body contributions to (Svjj) have the same general 
structure as the corresponding contributions to (5v^ E ). The direct term and the interacting 
exchange term have been evaluated and are presented below. 



W s LC (K,dir) = 

£ sS^ / ^ { rfP ^- fP ') KPmp ' AP ' £ / Snk if!k)2 A9 (D " + Dm * + D ^ (3 ' 21) 



pmp'q 



W s LC (K,ij-ex) = ]T 



pmp'qm'n 



1 P 



8m 2 s 



(An) 171 \ 1 

dr ) ( .2' ' 



P / ^ 3 r lfe U7k? A q (D qm , + + X -D qn 



(3.22) 



The separable three-body F-diagrams have a more complicated structure for the 5vff. 
The direct diagram has the general form 



2 

d3r " d3r ' t 

pmp'qq' 



{/gOf,, /« 0» t } (r„ ■ V fi ) (V (vJOj) ■ V fl ) {/£og, /M}] . (3.23) 



The integrand can be written as a sum of four terms having the gradients in Sv^ c acting on 
different parts of the correlations: 



dv 1 



-\ ( fpl rT fP ) {f q ^ lk c [{°hoi}oT 3 {ot 3 Mr v)r y • vo^}] 



_ (/VT*^. [{Of,, 0« fc } (r« • V R )V R • (/JvO^ {Og, Qj})] (3.24) 
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The C-parts in the above expressions depend on the cosine of the azimuthal angle, 6i, and 
therefore cannot be expressed exclusively in terms of the D pq matrices, which include an 
implicit average over that angle. The C-parts must then be individually evaluated for each 
operator product. 

The contribution of the direct diagram is dominated by the first term in the above sum. 
In evaluating the interacting exchange and passive exchange diagrams, only the term corre- 
sponding to that dominant part has been included. For example, the interacting exchange 
contribution is approximated by 

pmp'qq'n \ / 11 



fik., 

r 



(/*)" - \ {f 9 ')' cos 2 6< + l - (/*)') \C [Ol {0% 01} O- {Og, Ot }] . (3.25) 



The corresponding passive exchange contribution was calculated with the additional sim- 
plification of considering only leading term contributions, having at least one central ff k 
link. 

We have also evaluated central-chain diagram contributions, denoted by W cc h, to the 
boost interaction expectation values. These diagrams are obtained by dressing two-body 
cluster diagrams with hypernetted central chains. Their contribution was found to be sig- 
nificantly smaller than the contribution from separable three-body terms. 

The combined results for the two-body clusters and for the many-body clusters appear 



in Table [VIII] . The first four columns contain the boost contributions to the energy of SNM 
and PNM, calculated using the optimal wavefunctions for the A18 interaction alone, while 
the next four columns contain the corresponding boost contributions, calculated with the 
optimal wavefunctions for the A18+UIX. As expected, the many-body contributions to the 
boost interaction energy are comparable to the two-body contributions, because Tp and 
{T)2B have similar magnitudes. 

Detailed breakdowns of the boost contributions in the A18+UIX* model are presented in 
Tables IX and |X| for 5v RE and Sv LC . In all cases the bulk of the contributions come from the 
direct terms. Since we can integrate (6v LC ) by parts to obtain (5v RE ) plus additional terms, 
their expectation values should be similar in magnitude. The ratio of the contributions of 
5v LC and 5v RE in SNM is found to be ~ 0.7 at all densities, while in PNM it is ~ 0.75 at 
p < 0.2 fm~ 3 , and ~ 1 at higher densities. Results of VMC calculations |37|] have shown 
that this ratio is ~ 0.5 in 3 H and 4 He. 

The Sv RE has two factors, P 2 /8m 2 and v s , suggesting the approximation 

= -L {T) {v,) - (3 - 26) 

Our results show that this approximation has errors of only ~ 10 % in SNM, but can be 
wrong by ~ 50 % in PNM. Nevertheless it may be used to estimate the order of magnitude 
of the contribution of relativistic boost interactions. 
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B. Nucleon Matter Energies 



We evaluate the boost interaction contributions as a first order perturbation. Thus the 
energies of nucleon matter with the A18+Sv interactions are obtained by simply adding the 
5v contributions listed in Table [VIII| to the A18 energies. The results are listed in Tables |VI 
and [VI % and shown in Figs. ^| and ||. 

The three-body interaction, UIX*, to be used with A18+5v, contains the term that 
is 0.63 times the V& in the UIX. Since the boost effects are treated only in first order, the 



energies for the Al8+5w+UIX* interaction are obtained by adding {[Svij — 0.37 VJL]) to the 



%1 

'%jk\ 

energies for A18+UIX interactions. The results are listed in Tables |VT| and [VTT| , and shown in 
Figs. § and |3|. At low densities the A18+UIX and Al8+5v +UIX* interactions yield similar 
results. However the energies predicted by the latter model are lower than those of the 
former at higher densities, where 0.37 (V^) is much larger than (5v). The difference between 
the energies predicted by A18 + 5v and Al8+5v +UIX* interactions at higher densities is 
due to three-body forces; it is smaller than that between the energies obtained from A18 
and A18+UIX by almost a factor of two. 

The results obtained with the Urbana density- dependent interaction (U14-DDI) [M are 
also shown in Figs. ^| and [3] for comparison. Since Skyrme-type interactions based on the 
U14-DDI E(p) explain nuclear binding energies quite accurately, it is likely to provide a 
reliable representation of phenomena at p < p$. 

The variation of the kinetic and interaction energies with nucleon density in our most 
realistic model, with Al8+5v+UIX*, interactions is shown in Figs. [| and |5|. Due to a large 
cancellation between the contributions of V 2n and V R *, the total Vy* contribution is now 
smaller than that of v^, and the 5v contribution is also small. An order of magnitude 
estimate of the boost correction to the three body interaction is given by (T)(Vijfc)/3m, 
using generalizations of Eq. ( |3.26| ). This estimate is less than 10 MeV at the highest density 
considered. 

The most significant remaining problem appears to be our neglect of the boost corrections 
to the v MD , the unboosted contribution of which is quite large at high densities. Such 
corrections involve terms of order (velocity) 4 , which are beyond the scope of the present work. 
The kinetic energy also has corrections of that order, which we cannot include because the 
A18 interaction is fitted to data using nonrelativistic kinetic energy. However, the correction 
to the Fermi gas kinetic energy is 3kp/56m 3 , which is less than 10 MeV in PNM at the density 
of 6p - 

The Dirac-Brueckner approximation provides another way to estimate relativistic effects 
using realistic models of nuclear forces. The energies of our Al8+5v +UIX* model are lower 
than those of lowest order Dirac-Bruekner calculations with the Bonn A potential [^TJ. For 



example, at 4po we obtain 58 and 128 MeV per nucleon, while the Dirac-Bruekner calculation 
gives 76 and 164 MeV per nucleon, for SNM and PNM respectively. However, the difference 
between the results of the two methods is significantly smaller than that between either of 
the two and the uncorrected A18 energies of 3 and 59 MeV per nucleon. (Tables |VT| and 

ra- 



in the last column of Table |VJ we list SNM energies obtained by adding a correction 
72p 2 e _73P to the Al8+5v +UIX* results. The empirical binding energy and density of SNM 
are reproduced with 72 = 2822 MeV fm 6 , and 73 = 18.34 fm 3 . This correction has a 
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maximum value of 4.5 MeV at 0.11 fm -3 , and the "corrected" E(p) of SNM is a better 
representation of known nuclear properties at lower densities. It becomes identical to the 
E(p) obtained using the Al8+6v +UIX* model at higher densities. 



IV. COLD CATALYZED NUCLEON MATTER 

In this section we use the results of the earlier sections to calculate the equation of state 
and composition of cold, catalyzed matter, i.e. matter at zero temperature in its lowest 
energy state. Here the matter is assumed to be made up of neutrons, protons and leptons; 
the possible admixture of quark matter is considered later. Since the boost interaction is 
clearly an integral part of the two-nucleon interaction, we regard the models Al8+Sv and 
A18+5v +UIX* as realistic, and discuss their results in detail. The difference between the 
two models demonstrates the effect of the three-nucleon interaction. Some of the results 
obtained with the less realistic models without the boost interaction are also presented for 
comparison. 

Matter at zero pressure, at the surface of a neutron star, is made up of atoms of 56 Fe, just 
as in terrestrial iron. This is the most stable form of electrically neutral matter composed 
of neutrons, protons and electrons. Below the neutron star surface, as a result of the 
increased pressure of the matter caused by the gravitational attraction, the atoms become 
completely ionized and the electrons form a relativistic Fermi gas, whose Fermi energy 
becomes competitive in magnitude with nuclear energies ||42|| . Consequently, electron capture 
by the protons can occur. As a result, with increasing depth below the stellar surface the 
nuclei become more neutron rich, and cross the neutron-drip line. At this point the most 
energetic neutron orbitals have become unbound, and the matter consists of neutron-rich 
nuclei immersed in a neutron gas, whose density also increases as the pressure is increased 
further. The baryon number density at which this transition occurs is about 2 x 10~ 4 fm -3 . 
As the pressure and density continue to increase, the charge number Z of the nuclei remains 
in the range 30 — 40 but the mass number A grows steadily, and the distance between nuclei 
decreases. In the density range above about 0.06 fm~ 3 , where the volumes occupied by 
nuclei and by the surrounding neutron gas become comparable, the matter may undergo 
an inversion to bubbles of neutron gas surrounded by nuclear matter by going through a 
progression of phases involving non-spherical shapes, although this aspect of the crustal 
structure is somewhat model dependent [jTl) , [l3[| . At a density above about 0.1 fm -3 there no 
longer occur nuclei or other clumps of proton-containing matter, and cold, catalyzed matter 
becomes a uniform fluid of neutrons with a small fraction of protons. 

A reliable discussion of the properties of matter over the crustal density range requires 
a nuclear model that can describe inhomogeneous matter in the geometries occurring there. 
The difficulties encountered cause the model dependence mentioned. However, for neutron 
stars with masses £1.4M Q , the mass fraction contained in the crust of the star is less than 



about 2%. We have therefore used results of earlier work |TT| , |43[| for matter at densities 
< 0.1 fm~ 3 . Since over that density range our present matter energies agree well with ones 
used earlier, this substitution causes negligible inconsistency in our conclusions about the 
total mass of the neutron star. 

At densities of 0.1 fm -3 and greater, we require properties of charge neutral uniform 
matter made up of neutrons, protons, electrons and muons in beta equilibrium. At a given 
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baryon number density p = p n + p p , the conditions to be satisfied by the components are 
charge neutrality, 



and beta equilibrium, 



pp = pe+P^ (4.1) 



l^n pp Pei pfi pe- (4-2) 



Here /?i is the number density of the species indicated by the subscript i, and pi is its chemical 
potential including rest mass. After achieving a solution to these conditions, one obtains 
the total pressure, 



and total energy density 



P = P N + P e + P M , (4.3) 



e = e N + e e + Bp, (4.4) 



of the matter, and thereby the EOS, in the form P = P(p), where p = e/c 2 is the mass- 
energy density of the matter. The composition and other thermodynamic variables are also 
determined in this process. 

The baryonic part of the matter consists of strongly interacting neutrons and protons 
with a given proton fraction x p = p p / p. The minimum processing necessary in order to 
obtain an EOS for any of the sets of SNM and PNM energies of earlier sections is to fit the 
energies to a smooth function of density, so that the derivatives needed to obtain chemical 
potentials (/ij = de/dpi) and pressure (P = J2% PiPi — s) can be calculated. One then needs 
to interpolate between the x p = 0.5 SNM and x p = PNM results in order to find the value 
of x p required by beta equilibrium. 

The fitting procedure used earlier E3, and which we employ here, introduces more than 



the minimum processing just described. It uses a generalized Skyrme interaction containing 
a momentum and density-dependent delta function interaction. By a well-known procedure 



45], for an uniform system of infinite extent the matrix elements of this interaction can be 



written as an energy- density functional in the form 

H eff = (7^ + f(p, t p + (J^- + f(p, 1 - Xp)^ r n + g(p, Xp) + H gradient . (4.5) 

The number and kinetic densities p n (p) and t u m are defined in terms of the neutron or proton 
orbitals (f) a , n (p)(r) by pi(r) = J2a ni(a, T)\<p a ^(r)\ 2 and r»(r) = £ Q rn(a, T)|V0 Qji (r)| 2 , where 
the Fermi-statistics density of states rii(a, T) = [1 +exp((e Q . ) i — p^/T)]^ 1 at T = become a 
step function 0(e a) j — p^. The function / in Eq. (|4.5|) comes from the assumed momentum 
dependence of the interaction, and is parametrized in the form 

f(p,x P ) = f (p) +x p f 1 (p). (4.6) 
It produces density- and isospin- dependent effective nucleon masses m*(p,x p ), 
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+ f(p,x p ), ——, - = — + f(p,l-x p ). (4.7) 



2m*(p, x p ) 2m ' ' 2m* (p, x p ) 2m 

The term H gradient in Eq. (|4.5|) depends on gradients of the densities, and thus affects finite 
systems; it is determined uniquely by the functions / f43|. The relation between n and pi 



depends on temperature via the Fermi functions. Therefore knowledge of the matter energies 
at finite temperatures, as was given in Ref. flCl , permits a determination of f(p,x p = .5) 



and f(p,x p = 0), i.e. fo(p) and fi(p) in Eq. (|4.6j). The potential-energy terms g(p } x p = .5) 
and g(p,x p = 0) may then be obtained by comparison with the zero-temperature energies. 
The functional forms used for f(p, x p = .5, 0) and g(p, x p = .5, 0) are chosen to represent 
with appropriate accuracy the calculated energies; it is not sufficient, nor necessary, that 
they have the extremely simple power dependence on p of customary Skyrme interactions. 
The effective interaction FPS obtained in this manner from the SNM and PNM energies 



calculated with the U14-DDI interaction has been described and used elsewhere |4^j43|JTT 

In the present work, we have studied the matter energies only at zero temperature, and 
cannot therefore make a new determination of the momentum-dependent (r-containing) 
terms in H e ff. Rather than omitting them, we have kept intact the r-dependent terms from 



the earlier determination and have modified only the potential energy term g(p, x p ). To 
facilitate handling of the two different phases, we make separate fits to the normal low density 
phase (LDP) and the high density phase (HDP) with pion condensation. The analytic forms 
used for the fitting are given in Appendix A. They are chosen solely to provide an economical 
fit to the calculated energies. 

The interpolation between x p = 0.5 and x p = is carried out assuming a (1 — 2x p ) 2 
dependence of the energy at a given density. It is well known that matter energies, as a 
function of be expanded in powers of (1 — 2x p ) about x p = 0.5. Previous studies 

4l3fl , using cluster expansion techniques 0, have found that the quartic terms are small, and 



that the quadratic terms determined from results at x p = and 0.5 are sufficient to obtain 
few percent accuracy in the interpolation. The potential energy terms g in the effective 
interaction Eq. ( |4.5|) obtained by fitting to PNM and SNM are therefore interpolated by 
means of the expression 

g(p,x p ) = g(p,x p = .5) (l - (1 - 2x p f) +g(p,x p = 0)(1 - 2x p f. (4.8) 



A. The Phase Transition 

In fitting the results of the calculations described in Sections I- 1 1 1 , we find that the SNM 
and PNM energies coming from the models that include a three-nucleon interaction have a 
clear discontinuity in slope, associated with the phase transition, so that different analytic 
forms are needed above and below the critical densities, p t = 0.20 fm -3 for PNM, 0.32 fm -3 
for SNM. In Fig. |6| the energies and the fits are shown for the Al8+<5t> +UIX* model. (Here, 
and for the rest of the paper, we use for the SNM energies of this model the "corrected" 
values given in the last column of Table 0.) The Al8+5v and A18 models do not show 
such a discontinuity in any marked way. 

An important assumption of our treatment of the two phases of nuclear matter and the 
phase transition, as exhibited in Fig. || is that the analytic forms fitted to the PNM and 
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SNM energies of each phase may be extrapolated beyond the density region in which they 
are determined. Since the equilibrium is a two-dimensional phenomenon, the energy as a 
function of p and x p can be represented by a surface above the p, x p plane. The interpolation 
described in the last subsection makes the energy surface of the LDP a valley parabolic in 
the x p dimension extending from x p — to x p — 1, with its minimum at x p = 0.5. A similar 
interpolation between the energies of the HDP produces another parabola. Note that because 
of the charge symmetry breaking terms in the A18 interaction only the x p < 0.5 side of the 
parabolas is useful, and that the parabolic shape approximation is tested only in the LDP 

' I- 

A transition to a neutral pion condensed phase of PNM was obtained earlier with the 



older Argonne v 14 (A14) NNI and Urbana VII (UVII) TNI interactions by WFF fl|]. How- 
ever, with those interactions SNM has a normal ground state. WFF estimated the properties 
of the A14+UVII model of cold catalyzed matter by interpolating between the normal SNM 
and the pion condensed PNM. For the Al8+5v +UIX* interaction, we find that such a pro- 
cedure overestimates the symmetry energy of the HDP by ~ 10 to 20% in the density range 
0.2 to 0.5 fm~ 3 . Fortunately we see the transition for both PNM and SNM, and can avoid 
that problem. 

In Fig. [7] we show for the A18+5v +UIX* model the density at which the LDP and 
HDP E(p,x p ) surfaces intersect, i.e. where the interpolated phases have the same energy. 
The curve obtained is not necessarily a parabola in x p , since the density dependence of the 
curvature of the parabolas is more complicated than quadratic. Also shown is the proton 
fraction of beta-stable matter in each phase. 

The transition of matter from the LDP to the HDP probably occurs via a mixed phase 
region in which each phase need not be charge neutral as discussed in the following 
section. Here, for simplicity, we assume that each phase is a charge-neutral fluid and make 
a Maxwell construction to obtain the change in density due to the phase transition. The 
baryon chemical potentials, which for beta equilibrium are equal to the neutron chemical 
potentials, and total pressures of the two phases are equated to obtain the densities and 
proton fractions of the LDP and HDP in equilibrium. They are (0.204 fm~ 3 , 0.073) and 
(0.237 fm" 3 , 0.057) respectively. In the region between these two densities the most stable 
form of the matter is a mixture of the two phases. 



B. Proton fraction and \i e 



The proton fraction and the related electron chemical potential of the matter are impor- 
tant in assessing the cooling rates of neutron stars WJ\ and the possibility of kaon conden- 
sation in neutron star interiors j48"l , f49f . For the models discussed here we plot in Fig. |8| the 
proton fraction x p against baryon density. This figure also shows, for comparison, results 
obtained for other models and by other methods: the U14+UVII model, a predecessor of the 
present A18+UIX model, as given by WFF with the VCS method; the A18 model, results 



of Engvik et al. |T{| with lowest order Brueckner (LOB) calculations; the U14-DDI model 
|T0| using the FPS effective interaction (Tl|. Indicated on the plot is the critical value of x p , 



as a function of p, for the onset of the direct Urea cooling process in the presence of both 
electrons and muons. That process allows energy to be conducted from the interior of the 
star by neutrinos generated in binary thermal collisions, a very efficient process WTj . In the 
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models containing three-nucleon interactions, the discontinuity in x p at p ~ 0.2 fm~ 3 signals 
the onset of the HDP. The corresponding electron chemical potentials are plotted in Fig. || 

As suggested in Ref. ||12|| , the determining quantity for x p is the symmetry energy, given 
approximately by E sym ps E pnm — E S nm- The LOB and VCS methods yield similar values 
for x p with the A18 interaction up to p = 0.6 fm -3 . Beyond this density the x p obtained from 
VCS calculations starts to decrease, while that from LOB calculations continues to increase 
with density. The VCS and LOB energies for the A18 model are compared in Fig. |TU[ They 
are not too different for PNM. For SNM, however, at p > 0.6 fm -3 the LOB energy is much 
lower than that of the VCS. The VCS E sym saturates for the A18 model at p ~ 0.6 fm -3 , 
while that of LOB calculations continues to increase with density, causing the x p to do the 
same. It will be interesting to see if three hole-line and higher terms PU5l| neglected in the 
LOB calculations reduce the difference between LOB and VCS energies. 

Engvik et al. [[HJ have calculated x p in cold matter using the LOB method with all five 
modern NN potentials that provide high precision fits to the Nijmegen NN scattering data 
base. At a density of one nucleon per fm 3 the CD-Bonn model Jl7 gives the largest x p of 



0.15, while the Nijmegen I model [|T5] gives the smallest value of 0.10. The spread in these 
values is comparable with the difference between VCS and LOB results for A18 of x p = 0.09 
and 0.14 at this density. 

The 5v term and the three-nucleon interaction increase the symmetry energy, and push 
the x p barely above the Urea limit at high densities. For the Al8+5v +UIX* model the 
threshold is at a density of p = 0.78 fm -3 , and, as discussed in the next section, stars must 
have a mass >2.0M© to achieve such a density. However, this density is at the limit of 
our calculations and of the input physics. For example, admixtures of quark matter with 
hadronic matter, considered in the next section, may affect the Urea process in matter at 
such densities. 

The U14-DDI (FPS) model predicts values for x p that are much smaller than those 
predicted by all other models considered here, and in fact go to zero for p ~ 1 fm -3 . It 
is based on the U14 NN interaction, also used in the U14+UVII model. However, instead 
of adding the UVII three nucleon interaction to obtain the empirical saturation density 
of nuclear matter, it uses a density dependent modification (U14-DDI) of the U14 NN 
interaction |§ chosen to reproduce the energy, density and compressibility of equilibrium 
nuclear matter. Unlike the UVII interaction, this modification reduces the symmetry energy, 
and thus the x p , at high density. The main advantage of using three-nucleon interactions, 
instead of density dependent modifications of the two-nucleon interaction, is that the former 
can be tested via accurate calculations of the light nuclei. Unfortunately, the available 
results [3D] indicate that the UIX model may be overestimating the repulsion between three 
neutrons, thus overestimating the x p \ an improved version of the UIX model is currently 
being developed. 



V. NEUTRON STARS 

Using the methods just described we obtain for each model the EOS for cold, catalyzed 
beta-stable matter. At a baryon number density of 0.1 fm -3 they are joined onto an ear- 
lier EOS in which properties of the crust material has been treated more accurately p"T |. 
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The Oppenheimer-Volkoff general relativistic equations for a spherically symmetric (non- 
rotating) neutron star [|l| are 



dP (p + P/c 2 )G (m(r) + Anr^P/c 2 ) Mr) 

~dr = P ' mW 



rr 

/ Airr 2 p dr, (5.1) 
Jo 



where A(r) = (1 — 2Gm(r)/rc 2 )~ 1 . The corresponding equations for obtaining the moment 
of inertia, for a slowly rotating star, are given in Appendix B. Starting from some central 
mass-energy density p c , or equivalently from a central number density p c , these equations 
are integrated outwards to a radius r = R, at which P is zero, thus yielding the stellar 
radius, R, the gravitational mass of the star, m(R) = M, and the moment of inertia /. 
The dependence of the neutron star mass on central baryon density p c for the four models 



is shown in Fig. [11]. In order to estimate the effect of beta-stability on these results, we show 
also the trajectories obtained by using the pure neutron matter EOS for densities greater 
than 0.1 fm -3 , joined to the crust results of Ref. JTI|. Earlier results with the FPS EOS JTTJ 



are included for comparison. For the same set of results, the neutron star mass is plotted 



against the star radius in Fig. 12. 



The maximum masses for the five models illustrated in Figs. [IT] and [12] are listed in 



Table [XI]. While the models based on only two-nucleon interactions have maximum masses 



at or below 1.8M , those for the two models containing three-nucleon interactions have 
maximum masses well above 2M . The model that we believe includes most of the necessary 
physics is Al8+Sv +UIX*, which yields a maximum mass of 2.2M . This model achieves its 
maximum mass for a central baryon density p c = 1.14 fm -3 , which is not far beyond the 
calculated VCS energies, so that only slight numerical extrapolation is involved. 



The moments of inertia given in Table |XJ are the maximum values for each model. They 
occur for slightly lower central densities than do the maximum masses. The effect of the 
three-nucleon interactions on the maximum moment of inertia is in general considerably 
greater than on the maximum masses, because TNI tend to increase stellar radii also. 

The relativistic correction to the Fermi energies, — 3fc^/56m 3 , mentioned at the end of 
section III-B, when inserted (with its attendant contributions to chemical potentials and 
pressure) into the EOS of the model Al8+<5t> +UIX*, produces at p = 0.1 fm -3 a reduction 
in pressure of 6.9%, but at p — 1.0 fm -3 the reduction is only 1.1%. The effect on the 



maximum mass is a 0.2% reduction, smaller than the digits quoted in Table XI . 

The mass limits obtained with PNM EOS are within ~ 0.5% of the values determined 
from the EOS of matter in beta equilibrium, and the radii are within ~ 5%. Thus the 
uncertainties in the proton fraction of matter discussed in the last section do not have a 
large effect on these aspects of neutron star structure. 



A. Superluminality and Maximally Incompressible Matter 

Indicated on the curves in Fig. |TT] are the densities p s i at which the sound speed c s = 
J dPj dp becomes greater than the speed of light, c. Superluminal behavior would not occur 
with a fully relativistic theory, and it is necessary to gauge the magnitude of the effect it 
introduces at the higher densities. Kalogera and Baym provide one method for doing 
this. Following Rhoades and Ruffini ||53|| , they assume that the stiffest physically allowable 
EOS produces matter with a sound speed equal to the speed of light, i.e. 
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Since the matter is at zero temperature, for a single phase the partial derivative becomes a 
total derivative. Thus for higher densities than some pi the EOS is replaced by 

P(p) = P(p I ) + (p-p I )c 2 , (5.3) 

and the matter interior to pi is then maximally incompressible. Other thermodynamically 
related quantities may be obtained easily ||26||; the corresponding baryon number density is 



P = Pi 



P + UPi-Pic 2 ) (54) 



\ pj + \(P I -p lC ^ 



The relatively small reduction in the energy of the Al8+5i> +UIX* model resulting from this 
replacement for pj = p s i = 0.86 fm~ 3 is shown in Fig. 13. In Fig. 14 we show the effect on 



the mass vs radius plot. The maximum mass is reduced from 2.2OM to 2.19M , a very 
small change. In all models of the present work superluminal behavior occurs in stars only 
very close to the maximum mass limit, as can be seen from Fig. [IT]. Therefore replacing 
superluminal matter with maximally incompressible matter has little effect on the stellar 
properties. 

It is also possible that due to neglect of four-body and higher forces, and relativistic cor- 
rections of order (velocity) 4 and higher, the present work underestimates the sound velocity 
at lower densities. The effects of this possibility can be studied by assuming that the EOS of 
the Al8+<5t> +UIX* model is valid up to a chosen density pj, beyond which it is maximally 
incompressible. The results obtained for pi = 2, 3 and 4p are shown in Figs. and [TJ], 
and Table |Xj. The difference between the E(p) of the matter obtained by assuming that 
pi = 2po and the E(p) of the Al8+5v+UIX* model is several times the contribution of the 



UIX* interaction (see Fig. |14|). It therefore appears unrealistic to assume that pi can be as 
small as 2p . Results obtained with pi = 3p provide a better indication of what we can 
expect from the hardest EOS consistent with realistic models of nuclear forces. 



B. Transition to Quark Matter 

It is also possible that the present EOS is too hard due to the assumption that neutron 
star matter contains only nucleons and leptons. Should it also contain other hyperons such 
as A, £~ ,0 ' + and A _ ' 0,+ ' ++ , the EOS may be softer than that obtained with nucleons only 
0| . The forces between hyperons and nucleons and between hyperons are not as well known 
as nuclear forces, and it is therefore difficult to estimate whether such exotic species are 
present in neutron star matter. 

The chemical equilibrium in matter containing nucleons, lambdas, sigmas, deltas and 
leptons is governed by the equations 

Py,- = Pa- = Pn + Pe, (5.5) 
Pa = Pt,o = Pao = p n , (5.6) 

PY.+ = pA+ = Pp = Pn~ Pe, (5.7) 
PA++ = Pn~ 2Pe- (5.8) 
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If we neglect the interaction between nucleons and these hyperons, then the chemical po- 
tential of a hyperon at threshold density is given by its mass. Under such an assumption, 
the negatively charged hyperons appear in the ground state of dense matter when p n + p e 
reaches their mass, while the neutral hyperons appear when p n equals their mass. Since 
He in dense matter is larger than the mass differences between lambdas, sigmas and deltas, 
the S~ and A~ will appear at lower densities than the neutral A, provided the interaction 
effects are small. The chemical potentials of electons, and of neutrons in beta-stable matter, 



are shown in Fig. [L^ for the Al8+5f+UIX* and Al8+Sv models. The threshold densities 
for non- interacting A~ and A are indicated by horizontal line segments. Given their 
relatively low values, it is clear that in the absence of interactions these particles would 
be present in most neutron stars. Results obtained using crude models of the interactions 
between hyperons and nucleons, and between hyperons [f|] indicate that the S~ and A~ 
have the largest effect of all hyperons on the EOS; however the magnitude of the effect is 
very sensitive to the interaction model. 

The available hyperon- nucleon scattering data has been reviewed recently by de Swart, 
Maessen and Rijken pl|,|55] along with the status of one-boson exchange models of the 
interactions between hyperons. Additional information on A-nucleon interactions can be 
obtained from the measured A-nucleus binding energies. These indicate the presence of 



ANN three-body forces that are as strong as the three- nucleon interaction [56,57]. In view 



of these uncertainties, particularly concerning the important S~ and A~ interactions, we do 
not attempt to estimate the effect of these hyperons on the EOS of neutron star matter. 

A transition from hadronic to quark matter is expected at high densities. Knowledge of 
the EOS of both hadronic and quark matter is necessary to estimate the possible effects of 
this transition on neutron stars. Here, we use the present models of the EOS of hadronic 
matter, containing only nucleons and leptons, and the quark bag model with u, d and s 
quarks for the quark matter. The u and d quarks are taken to be massless, and s quarks 
to have a mass of 150 MeV. A Fermi gas of quarks of flavor i has density pi = kp^/ir 2 , 
due to the three color states. There is no one-gluon exchange interaction energy between 
quarks of different flavor, while that between quarks of flavour i is given by (2a/37r)Ei per 



quark % ||58|| . Here is the average kinetic energy per quark, and a is the strong interaction 
coupling constant, assumed to have a value of 0.5. The value of the bag constant B is poorly 
known, and we present results using two representative values of B = 122 MeV p^] and 
B = 200 MeV fm~ 3 @. 



The beta equilibrium conditions for charge neutral quark matter are 

+ He = = Hs, = /V (5.9) 

The energy densities of charge neutral quark matter and nuclear matter are plotted in 
Fig. fnjL In the interesting region of p ~ 1 fm~ 3 the total energy density of quark matter is 



about 1200 MeV fm~ 3 , of which only 122 or 200 MeV fm -3 comes from the bag. 

If, during the phase transition from nuclear to quark matter, the nuclear and the quark 
phase are each required to be charge-neutral beta-stable fluids whose pressures and baryon 
chemical potentials are equilibrated, then for the Al8+5w+UIX* model the transition is 
found to extend over the density range p = 0.86 <-> 1.57 fm -3 for the B = 200 MeV fm -3 
case. For B = 122 MeV fm -3 , the range is 0.79 <-> 1.20 fm -3 . The matter within this 
density range, as so treated, is a constant-pressure mixture of that at the two ends of the 
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range. Such a constant-pressure mixed phase does not occur in the neutron star: as can be 
seen from Eq. ( |5.1|) , since the pressure does not change, the density changes discontinuously 
from the lower nuclear matter density to the higher quark matter density. 

Since for the nuclear component of the matter at p ~ 1 fm -3 the electron screening 
length is ~ 7 fm, while the quark matter component has negligible electron distribution, 
the previous assumption that each fluid retains its electron distribution during the phase- 
mixing is not a good approximation at length scales ^7 fm. As described in Ref. |2"3"| , |2T)|] , 
an alternative assumption that the nuclear and the quark matter share a common electron 
distribution in the mixture is more appropriate. For a given nucleon matter density p^ the 
transition now involves equilibrating the chemical potentials, 

Pn = Pu + 2Pd, Pp = 2Pu + Pd, (5.10) 

and the hadronic pressures of nuclear and quark matter. In the process, the proton fraction 
of the nucleon matter and the densities of the quarks are also determined. The lepton 
chemical potentials are then given by 

Pe = Pfi = Pn - Pp = Pd - Pu, (5-11) 

from which the electron and muon densities are readily found. The fractional volume, J'q, 
occupied by the quark component in the mixture is now chosen so that the nuclear, quark 
and lepton components are in sum charge neutral: 

Pe + Pn = ^fq^pu - Pd- Ps) + (1 - /q)p p , (5.12) 
after which the macroscopic baryon density of the mixture can be calculated: 

P = Pmix = ~/q(p« + Pd + Ps) + (1 - fo)pN- (5.13) 

The procedure is carried out provided that < fq < 1. Some numerical details are given 
in Appendix C. 

The mixture of A18+5t> +UIX* model nucleon and the B = 200 MeV fm~ 3 quark matter 
occurs over the density range p m i X = 0.74 <-> 1.80 fm~ 3 . Over this range, a charge neutral 
mixture of quark and nucleon matter is more stable than the constant-pressure mixture of 
separately neutral nuclear matter and quark matter. The quark matter fraction Jq, the 
charge per baryon of the nucleon matter, given by the proton fraction x p , the charge, x c , 
and strangeness, x s , per baryon of the quark matter and the charge xg per baryon carried 



by the leptons in the mixed matter are plotted for this case in Fig. 17. 

Over the density range occupied by the mixture the nucleon matter has densities p^ = 
0.74 «-> 0.96 fm -3 , while the quark matter densities are 1.12 <-> 1.82 baryons fm -3 . The 
mixture thus consists of dense negatively charged quark matter immersed in somewhat less 
dense positively charged nucleon matter. The neutralizing charge of the leptons decreases 
in magnitude from its value in the purely nucleon matter, and tends rapidly to zero by the 
time the quark fraction has reached 50%. 

The lower density end of the mixed phase region is more relevant for neutron stars. Here 
one expects small drops of dense quark matter with x c ~ x s ~ — 1 appropriate for matter 
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made up of X~ hyperons. Ground states of dense matter can have such a form also if matter 
made of S~ has a softer EOS than nucleon matter. Knowledge of the interactions between 
£~ hyperons is necessary to further explore this possibility In regard to the direct Urea 
cooling process, the proton fraction of the nucleon matter in the mixture increases with fg, 
but the lepton fraction decreases. The required momentum balance among neutrons, protons 
and electrons is only achieved at a baryon density of 0.86 fm -3 , somewhat larger than the 
value 0.78 fm -3 needed for the nucleon matter alone. At densities larger than 0.86 fm -3 the 
momentum balance remains possible, although the decreasing density of electrons reduces 
the rate of the Urea process. Thus the indirect effect of the quarks on the nucleons is to 
delay the start of the Urea process somewhat. At a density of 0.86 fm -3 the quarks, which 
are about 7% of the matter, can if behaving as free particles contribute to the Urea process 
directly. Such considerations, however, are complicated by effects of quark matter occurring 
as droplets. 

At the high density end of the mixed phase, the quark matter has approximately equal 
number of u, d and s quarks, and nearly symmetric nucleon matter occupies a small fraction 
of the total volume. Here one expects small drops of SNM, i. e. nuclei compressed to a density 
of ~ 6po by the pressure of surrounding quark matter. However, this interesting form of 
matter at very large densities does not seem to occur in stable neutron stars according to 
the present calculations. 



The two-component equilibrium suggested in Ref. [23] is formally the same as the equi 



librium between nuclear matter and dripped neutrons in crustal neutron-star matter. Ex- 
tending that similarity, Ref. [24j includes the effects on the phase transition of the energy 
of the surface dividing nucleon and quark matter, as well as the Coulomb energy associated 
with the difference between the charge densities of nucleon and quark matter. There then 
occur in the mixed region varying sizes and shapes of matter in one phase surrounded by 
liquid in the other phase. The presence of a new parameter of unknown magnitude, the 
surface energy, complicates the problem, however, and lacking new information on it, we 
have not included surface and Coulomb effects in results given here. 

The energy densities of the mixed region, calculated neglecting its surface and Coulomb 



contributions, are shown in Fig. 16. As discussed at length in Ref. [23|, the total pressure 



obtained is now not constant over the density range of the mixture, although as can be 
deduced from Fig. [16], it is reduced from that of the pure nuclear matter. The resulting 
neutron star masses are shown in Fig. [18] as a function of central baryon density. Maximum 
masses are given in Table DGI. For the Al8+5v model, admixing quark matter has small 
effects on neutron star structure; however, the effects are noticeable for models containing 
three-nucleon forces. Note that by neglecting the surface and Coulomb energies we underes- 
timate the energy density of mixed matter and thus overestimate the effect of the admixture. 
The maximum densities in stable neutron stars remain below those for pure quark matter 
in all the four cases considered. 



C. Transition to spin-ordered phase 

When each phase is required to be charge neutral, the mixed region between the normal 
nuclear matter LDP and the spin-ordered, or neutral pion condensed HDP has constant 
pressure and constant average nucleon chemical potential. In the neutron star such a mixed 
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region would not occur, and the density at this pressure would change discontinuously by 
0.033 fm~ 3 for the Al8+5v +UIX* model, a 15% jump. However, we should use the approach 
described in the last subsection for this transition also, and consider mixtures of charged 
LDP and HDP in a common lepton sea. On applying the equilibrium conditions of nuclear 
pressure and neutron and proton chemical potentials between the LDP and HDP, one finds 
that the density limits between which a mixture occurs is only slightly extended compared 
to those with the charge-neutral-fluid equilibrium quoted earlier, but the mixture does now 
have a pressure that changes with density. However, due to the similarity in charge character 
of the two phases the change in pressure in going from pure LDP to pure HDP is rather 
small, and therefore the thickness of the LDP-HDP mixed region in the star is also fairly 
small. 



Selected neutron star density profiles for the Al8+£t>+UIX* model are shown in Fig. IS . 
Although it is not detectable from the figure, the mixed region associated with pion con- 
densation now extends over a finite region of the star. Its radial thickness, Aty, is ~ 14 m 
for the 2.10M Q case, and ~ 40 m for the 1.41M case. The 2.10 and 2.00M Q stars shown 
in Fig. [19] have the same central density of 0.86 fm -3 ; their mass difference is due to the 
admixture of quark matter considered only in the 2.00 M star. The 1.41 M Q stars cannot 
have quark matter admixtures in the present models. 

It is interesting to know the structure of matter having a mixture of LDP and HDP. One 
possibility is that suggested for the nuclear matter-quark matter transition by Heiselberg, 



Pethick and Staubo [p4}| , that there occur in the mixed region varying sizes and shapes of 
matter in one phase surrounded by liquid in the other phase. If we apply their estimates 
to this phase transition, assuming a surface tension for the interface in the range a = 1 to 
10 MeV fm -2 , we find that the characteristic length scale a of the entities involved is ~ 60 
to 120 fm. The electron screening length, which for this approximation to be valid needs to 
be also of this magnitude, is only ~ 10 fm, however. This indicates that neither assumption 
about the mixture, that of two charge-neutral fluids (which results in Arj- = 0) nor of two 
nuclear components with a common electron fluid (which gives the Aty values given above) 
is correct, and we need to look at a mixture that is somewhere in between these extremes. 
That problem will be pursued elsewhere. 



VI. ADIABATIC INDEX T OF NEUTRON STAR MATTER 

A measure of the stiffness of matter described by the EOS P = P(p), is the adiabatic 
index T: 



p dP p o 



P dp P 



<■ (6.1) 



If T were constant, then the EOS would become P cx p r . In this form, it is called a polytrope, 
an idealized EOS on which many pioneering studies of stellar structure were based ||61|| . The 
values of T in several limiting cases are well known. At lower densities, the rest mass of the 
constituents dominates the energy density of matter, i.e. p ~ p m. In this limit, for matter 
without correlations, T = 5/3, 2 or 3 when the pressure is provided by Fermi kinetic energy, 
static two-body or three-body interactions respectively. Repulsive momentum dependent 
interactions lead to larger values of T. In the opposite extreme of very high density, the 
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energy per particle is much larger than the particle's rest mass. Neglecting rest mass, it 
is given by Xp n , where n = 1/3,1 or 2 for the above three cases, and A is a constant. 
Thus p = Xp n+1 , P = n\p n+l , and T = 1 for all values of n. This does not contradict 



the well-known value of T — > 4/3 in white dwarfs | 61fl . In white dwarf matter the nuclei 
are nonrelativistic and the electrons are relativistic. The energy density of this matter 
pc 2 = pMc 2 , where p is the number density of nuclei, and M is their mass. The pressure 
given by relativistic electrons is proportional to p 4 / 3 or equivalently p 4//3 , giving T = 4/3. 

Plots of r and the sound speed for various cases previously described are given in Fig. |20| . 
The two curves corresponding to pure nucleons behave as expected from the above limits. 
From the crust, T rises to ~ 3: of the two maxima, the greater value corresponds to the 
model with the explicit three-body interaction UIX*; T then decreases at higher densities. 
The small peaks followed by cusps in these curves at p/m ~ 0.15 fm~ 3 mark the softening 
of the EOS due to the opening of the muon channel. The deep discontinuity in both F and 
c s for the Al8+5u+UIX* model are due to the LDP-HDP transition. The indicated values 
T = 0.075 and c s /c = 0.045 shown at the bottom of the declevity are averages resulting from 
the treatment of that transition in Sec.V-C. (When treated in terms of a mixture of neutral 
fluids, in which case the pressure is constant, these quantities would be zero.) 

It is simple to express in terms of c 2 , the extent of the stellar radius Ar-p occupied by 
this phase transition from Eq. (157 



Arr „ Asa, ST „ (6.2) 
P9t r l T 

where p is an average quantity for the two phases, m(rT) is the mass interior to the radius tt 
at which the transition ocurs, A(r^) is the redshift factor there, and Sp ~ 8pm is the density 
range of the transition, given earlier. The gravitational acceleration at the transition 
region is greater for the more massive star, and explains the trend of our values obtained 
from the star profiles. 



The use of maximally incompressible matter beyond p = p s i, occurring in Fig. |20] at 
p s i/m = 1.06 fm -3 , produces a sharp reduction in T. This modification produces so little 
effect on the stellar mass limits (shown in Fig. |TB|) because it occurs at a density close to 
the central density in the star with maximum mass. The curves corresponding to nucleon- 
quark mixed phases discussed in Sec.V-D reduce T and c s more gradually: for an ultra- 
relativistic quark gas c s = cj v3, and the sound velocities of models with quarks are seen 
to be approaching that value. Note that the partial derivative involved in the calculation 
of r and the sound velocity of nucleon-quark mixed phases must be carried out holding the 
quark fraction Jq constant. 



VII. CONCLUSIONS 

The results obtained with the modern A18 NNI for the gross properties of neutron stars, 
such as the mass limit and the radius of 1.4M stars, are not very different from those 



obtained from the earlier U14 and A14 models ||12|| . The main difference is that with A18 
and A14 NNI there is the indication of the neutral pion condensation phase transition, while 
there is none with U14. However, the fact that A18 provides an exact fit to the NN scattering 
data while U14 does not weighs against the U14 results. 
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As is well known, the cooling of neutron stars is accelerated by pion condensation E7 



The detailed mechanism associated with the spin-ordered condensate that we find is as yet 
unexplored. In addition, we estimate that due to this phase transition, there may be regions 
in the star where the matter density changes rapidly over a ten meter distance scale. The 
matter in this thin layer of the star can have interesting structures on a ten fermi length 
scale. 

We find that inclusion of the relativistic boost correction, 5v, to the NNI increases the 
mass limit from 1.67M to 1.8M without any TNI, while it is reduced from 2.38M Q to 
2.2OM upon inclusion of the Urbana models of TNI. The reduction occurs because the 
TNI needed to fit nuclear data have weaker repulsive parts after including 5v. For the same 
reason, the effect of Urbana models of TNI on the mass limit is diminished from 0.71M Q to 
0.4M Q when the Sv is included. Note that these models of TNI have only two terms; the 
long range two-pion exchange term, and a short range term with no spin-isospin dependence. 
Their strengths are determined from the density of SNM and the triton energy. However 
the triton, which has isospin 1/2, is insensitive to the interaction between three neutrons. 
Improved models of TNI must consider data such as the binding energy of 8 He, which are 
sensitive to the interaction between three neutrons. 

If the effective value of the bag constant B is larger than 122 MeV/fm 3 it appears that 
only the heaviest neutron stars may have small drops of quark matter in their interior. The 
quark composition of these drops is similar to that of aggregates of S~ hyperons. If the 
interaction of S~ with dense nucleon matter is not repulsive, there may exist S~ hyperons 
in nucleon matter at densities below the threshold for the appearence of quark drops. It is 
necessary to build AE~ and NNH~ interaction models to more fully explore this possibility. 

Using maximally incompressible matter at high densities we find that the upper limit 
for the maximum mass consistent with nuclear data is ~2.5M Q ; this is not far from the 
prediction of 2.21M Q of our Al8+5t> +UIX* model. The lower limit for the maximum mass is 
more difficult to establish due to the unknown interactions between hyperons and nucleons. 
It may be as low as 1.74 M® if the quark bag constant B has a value of ~ 122 MeV/fm 3 , 
and if the three neutron interaction is not as repulsive as the Urbana model IX at high 



densities, even without admixture of hyperons in matter. Recently several authors [|62]-|64 
have argued that there are indications of the existence of neutron stars with M ~ 2M Q . If 
these are confirmed, then models without TNI will be ruled out. However, such a possibility 
is still the subject of active debate [BH . 
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APPENDIX A: FORM OF EFFECTIVE HAMILTON! AN 



The effective interactions that fit the models examined in the first part of the paper all 
have the form 

H = ^ + (p 3 + (1 - x p )p 5 )pe- pip) j r n + (J^ + (p 3 + x pP5 )pe- p4p] j r p 

+ g{p, x p = .5) (l - (1 - 2x p f) + g{p, x p = 0)(1 - 2x p ) 2 , (Al) 
where p = p n + p p , and at zero temperature 

r p = 3 -(3n 2 p) 2 / 3 xf , r n = jj(3vr 2 p) 2/3 (l - x p f\ (A2) 

The parameters defining the r-dependent terms are the same for all of the models, and are 
given in the caption of Table [X1J| . For the A18 and Al8+<5t> models at all densities and for 
the LDP of models with TNI, the parametrization is 

g L (p, x p = 0.5) = -p 2 (pi + p 2 p + Pep 2 + (pxo + Pup) e" p s p2 ) , 
g L (p, x p = 0) = -p 2 (p 12 /p + p 7 + p 8 p + Pis e" pip2 ) , (A3) 

while for the HDP of models with TNI, 

9h(p, x P = 0.5) = g L (p, Xp = 0.5) - p 2 (pn(p - pig) 

+ P2i(p-Pi 9 ) 2 )e pi8( ^ pi9) , 
9h(p,x p = 0) = g L (p,x p = 0) - p 2 (p 15 (p - P20) 

+ Pu(p-P2o) 2 )e pi6{p - p2o) , (A4) 



The values of the parameters are as given in Table XII 



APPENDIX B: MOMENT OF INERTIA OF A SLOWLY ROTATING STAR 

We use the general relativistic equations for a slowly rotating star as described by Hartle 
|66| . The metric for the nonrotating star is 



ds 2 = -e v{r) dt 2 + e x{r) dr 2 + r 2 [dQ 2 + sin 2 Odcf) 2 ). (Bl) 

It involves the radial functions u(r) and A(r). The Oppenheimer-Volkoff equations for the 
pressure P(r), mass function m(r) and A(r) are given in Sec. V, and A(r) = Zn(A(r)). The 
function u(r) is defined by 

dv _ 2G(m + 4vrr 3 P/c 2 )A(r) 
dr r 2 

with the boundary condition e u ^ R > = 1/A(R), and there is also an equation for the rotational 
drag, Lj(r), 
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d ( 4 .dw\ 3 dj _ 

dr \ dr ) dr 



(B3) 



Here j(r) = e~^^ r ' + ^ r ''' 2 ; it has the boundary value j(R)= 1. In the limit of slow rotation, 
such that the angular velocity Q <C GM/R 2 c, u{r) has the boundary condition u(R)/Q = 
1 — 2GI/R 3 c 2 . I is the total moment of inertia, given by either of the integrals 

2c 2 f R 3 dj(r)cu(r) , 8tt f R 4 / P\ . . . x cu(r) , 

This set of equations, together with (|5JJ), is integrated from r = to the value r = R where 
the pressure becomes negligible, with a given equation of state P = P{p), and a central 
density p(0) chosen to give the desired neutron star mass. One then has also the radius and, 
after satisfying the boundary conditions, the moment of inertia. 



APPENDIX C: DETAILS OF THE NUCLEON-QUARK EQUILIBRIUM 

The nuclear matter-quark matter strong interaction equilibrium discussed in Sec.V-B is 



illustrated graphically in Fig. [21]. The x p > (right side) of the plot relates to nucleon 
matter and the left relates to quarks. Although states of quark matter can extend to x c = 2 
only the x c ^0 part is relevant here. The chemical potentials of isobars corresponding to ten 



equally spaced pressures are shown as dashed lines. Equilibria corresponding to Eq. (]5.10|) 



and to pressure are represented by the rectangles. The top (bottom) horizontal sides of the 
rectangles join states with the same \i n on quark and nucleon isobars corresponding to 
the same pressure. (Only in a few cases do the pressures involved correspond to those of the 
displayed grid.) The vertical sides ensure that the chemical potentials, \i n and /i p on nucleon 
side, and 2/i rf + \i u and 2/i„ + \i d on the quark side, belong to matter with the same charge 
per baryon. Other properties of interest, such as densities, are not represented on the plot. 
These equilibria are a property purely of the strong interactions, irrespective of leptons and 
charge neutrality. The latter property then determines the quark fraction /q, defined in 
Eq. (|5.12| ). The rectangles labelled b, c, d and e have values of /q going from 0.001 to 0.999 



and represent physical states, while those labelled a and / give /q = —0.169 and /q = 2.38 
corresponding to unphysical states. However, they all correspond to equilibrium under the 
strong interactions. 
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FIGURES 

FIG. 1. Central and tensor healing distances (d c and dt), for SNM and PNM, with A18 
interaction alone (upper graph), and with A18+UIX interaction (lower graph.) 



FIG. 2. The energy per nucleon, E(p), of SNM for various interaction models. 
FIG. 3. The energy per nucleon, E(p), of PNM for various interactions models. 
FIG. 4. Kinetic and interaction energies in A18+<5t>+ULX* model of SNM. 

FIG. 5. Kinetic and interaction energies in Al8+5v +UIX* model of PNM. 

FIG. 6. The PNM and SNM energies for the Al8+5t>+UIX* model, and the fits to them using 
an effective interaction. The full lines represent the stable phases, and the dotted lines are their 
extrapolations. 

FIG. 7. On a plot of proton fraction x p vs. baryon density, for the Al8+5v +UIX* model, the 
boundary between the LDP and HDP, obtained in the manner described in the text. The dashed 
curve is the proton fraction of beta-stable matter, and the dotted lines mark the boundary of the 
mixed phase region. 

FIG. 8. For beta-stable matter, the proton fractions x p for the four models discussed in 
the text vs. baryon density. The dashed curve, U14-DDI (FPS), is from [O, the dotted line, 



A18(LOB) is from jig], the points are from []12|,[50|] for the U14+UVII model, and the dotted line 



is the threshold for the direct Urea cooling process, as a function of p. 

FIG. 9. For beta-stable matter, the electron chemical potentials p e for the four models dis- 
cussed in the text vs. baryon density. The dashed curve, U14-DDI (FPS), is from pT| , and the 
horizontal line, at p e = m^c 2 , is the threshold for a muon contribution to the lepton fraction. 

FIG. 10. Comparison of the energies of PNM and SNM obtained for the A18 model with the 
VCS method in the present work, and with LOB calculations by Engvik et al. [pTsf] . 

FIG. 11. Neutron star gravitational mass, in solar masses, vs. central baryon density, for 
the four models described in the text. The full curves are for beta-stable matter, and the dotted 
lines are for pure neutron matter. The vertical lines show the density above which the matter is 
superluminal. The dashed curve, FPS, is from jll]]. 
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FIG. 12. Neutron star gravitational mass, in solar masses, vs. radius, in kilometers, for the 
four models described in the text. The full curves are for beta-stable matter, and the dotted ones 
are for pure neutron matter. The dashed curve, FPS, is from [11]. 



FIG. 13. For beta-stable matter according to the A18+5v-|-UIX* model, the total energy per 
baryon vs. baryon density (full curve). The dashed curves are for the assumption that matter is 
maximally incompressible for densities greater than the indicated value, and the lower full curve is 
for the Al8-\-5v model. 

FIG. 14. Neutron star gravitational mass, in solar masses, vs. radius for the Al8+fe+UIX* 
model (upper full curve), for the maximally incompressible modifications of this model at densities 
beyond choosen values of pj (dashed curves), and for the Al8+5v model (lower full curve). 

FIG. 15. The neutron and electron chemical potentials in beta stable matter according to 
models A18+5v+XJIX.* (full line) and Al8+5v (dashed line). Threshold densities for the appearance 
of noninteracting hyperons are marked by horizontal line segments. 

FIG. 16. For beta-stable matter, the energies per unit volume for the A18+<fo+UrX* and the 
A18+5v models, and the quark bag models with B = 122 and 200 MeV/fm 3 are shown by full and 
dashed lines; the dotted lines correspond to neutral mixtures of charged nuclear and quark matter. 

FIG. 17. For the A18+5v+\JlX* models and the quark bag model with B = 200 MeV/fm 3 , 
the proton fraction of nucleon matter x p , the charge (x c ) and strangeness (x s ) per baryon of quark 
matter, the volume fraction /q occupied by quark matter, and the charge xe per baryon carried 
by the leptons in the mixture. 

FIG. 18. Neutron star gravitational mass, in solar masses, vs. central baryon density, for the 
A18+5f+UIX* and A18+5v models with (dashed curves) and without (solid lines) quark matter 
admixture. The two dashed curves for the Al8+<5t> model are for B = 122 and 200 respectively. 

FIG. 19. Density profiles of 2.1M Q and 1.41M Q stars of beta-stable matter using the 
A18+5i>+UFX* model without quark matter admixture (solid lines), and of 2.OOM0 star with 
quark matter admixture (dashed line), assuming a bag constant B = 200 MeV fm~ 3 . 

FIG. 20. The sound speed and the adiabatic index T vs. p/m for beta-stable matter according 
to the Al8+5v+\JYK* (solid curves) and A18+8v (long dashed curves) models. The short dashed 
curves are for the assumption that matter is maximally incompressible for densities greater than 
the indicated value. The dotted curves are for the nucleon-quark mixed phase. 
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FIG. 21. Graphical representation of the strong interaction equilibrium between the A18+ 
(5i>+UIX* model of nucleon matter and the B = 200 MeV quark matter for the density region 
of the mixing (see text). Dashed curves labelled 1 through 10 are strong interaction isobars, i.e. 
without considering charge effects and leptons. They are for pressures from 100 MeV fm -3 to 
550 MeV fm -3 in intervals of 50 MeV fm~ 3 . The upper (lower) set of curves on the right side 
show the neutron (proton) chemical potential vs proton fraction (right side) , and the quark-matter 
counterparts on the left side show fi u + 2fj,d (2fj, u + /id) vs charge per baryon x c . The corners of the 
rectangles correspond to the equilibrium values of the chemical potentials and charge fractions. 
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TABLES 



TABLE I. Cluster contributions to A18 SNM E{p) in MeV. 



p 


T F 


(T)2B 


(T)mb 


AT 


(V")2B 


U,MD\ 


(v s )mb 


l„,MD\ 

\V /MB 


5E 2 b 


0.04 


8.77 


5.33 


-0.29 


0.16 


-19.46 


0.03 


1.43 


0.17 


-0.27 


0.08 


13.93 


10.30 


-1.06 


0.47 


-36.84 


0.43 


4.44 


0.64 


-0.57 


0.16 


22.11 


19.64 


-1.79 


1.11 


-66.48 


2.04 


8.44 


2.58 


-1.13 


0.24 


28.97 


26.95 


-2.47 


1.49 


-91.46 


4.69 


12.06 


5.39 


-1.73 


0.32 


35.09 


34.45 


-3.92 


1.63 


-115.82 


7.89 


17.17 


9.37 


-2.37 


0.40 


40.72 


40.55 


-3.23 


1.53 


-135.38 


11.96 


17.57 


14.48 


-3.03 


0.48 


45.99 


46.22 


-1.95 


1.17 


-152.83 


16.69 


16.55 


20.76 


-3.63 


0.56 


50.96 


51.95 


-1.82 


0.46 


-170.59 


21.72 


17.92 


28.24 


-4.17 


0.64 


55.71 


58.03 


-3.73 


-0.99 


-188.73 


26.61 


22.56 


36.91 


-4.60 


0.80 


64.64 


68.60 


4.47 


-4.07 


-202.42 


34.67 


8.42 


51.76 


-5.13 


0.96 


73.00 


79.71 


-0.40 


-8.27 


-236.38 


47.43 


18.47 


80.12 


-5.43 



TABLE II. Cluster contributions to A18 PNM E(p) in MeV. 



p 


T F 


(T)2B 


{T)mb 


AT 


(V S )2B 


\V )2B 


{v s )mb 


{V /MB 


5E 2B 


0.04 


13.92 


5.08 


-1.93 


0.14 


-15.79 


0.41 


4.90 


-0.14 


-0.39 


0.08 


22.09 


8.21 


-2.49 


0.26 


-26.87 


1.34 


6.92 


-0.22 


-0.46 


0.16 


35.07 


15.00 


-3.02 


0.43 


-48.10 


4.15 


9.75 


0.28 


-0.80 


0.24 


45.95 


21.67 


-2.89 


0.20 


-67.45 


7.67 


11.01 


1.92 


-1.20 


0.32 


55.67 


28.50 


-3.65 


-0.54 


-86.52 


11.95 


12.81 


5.16 


-1.73 


0.40 


64.60 


31.71 


-0.28 


-1.07 


-98.77 


17.57 


6.97 


9.95 


-2.33 


0.48 


72.95 


37.39 


4.57 


-1.23 


-113.39 


23.89 


0.93 


15.57 


-3.00 


0.56 


80.84 


43.86 


3.15 


7.25 


-141.33 


35.90 


-1.98 


32.50 


-3.86 


0.64 


88.37 


56.20 


0.99 


8.55 


-168.18 


40.75 


0.50 


45.53 


-4.80 


0.80 


102.54 


67.08 


2.35 


5.72 


-201.26 


57.99 


-6.50 


73.07 


-7.00 


0.96 


115.80 


80.44 


-2.75 


2.83 


-237.50 


74.78 


-8.05 


112.66 


-10.08 
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TABLE III. Cluster contributions to A18+UIX SNM E(p) in MeV. 



p 


T F 


( T )2B 


(T") 


AT 


(V S )2B 


l„,MD\ 
\V >9R 


{v s )mb 


l„,MD\ 

\V /MB 


5E/2B 


0.04 


8.77 


5.70 


-0.34 


0.19 


-20.00 


-0.01 


1.60 


0.19 


-0.28 


0.08 


13.93 


10.16 


-0.97 


0.46 


-36.51 


0.44 


4.21 


0.64 


-0.60 


0.12 


18.25 


15.27 


-0.74 


0.89 


-51.63 


1.24 


5.25 


1.60 


-1.10 


0.16 


22.11 


20.97 


-0.81 


1.21 


-66.74 


2.12 


6.52 


3.01 


-1.80 


0.20 


25.65 


25.84 


-0.46 


1.51 


-80.18 


3.47 


7.12 


4.77 


-2.54 


0.24 


28.97 


30.51 


0.22 


1.87 


-93.16 


5.38 


7.51 


7.10 


-3.33 


0.32 


35.09 


39.38 


1.32 


2.46 


-117.90 


9.41 


8.60 


12.56 


-5.21 


0.40 


40.72 


53.32 


0.54 


6.24 


-154.44 


15.18 


23.65 


25.65 


-6.93 


0.48 


45.99 


61.13 


1.08 


7.80 


-178.46 


18.79 


27.12 


34.34 


-7.93 


0.56 


50.96 


69.68 


2.27 


8.51 


-201.66 


24.75 


27.70 


45.69 


-8.67 


0.64 


55.71 


79.46 


1.92 


10.20 


-227.19 


28.60 


32.93 


58.95 


-9.20 


0.80 


64.64 


93.07 


1.71 


12.89 


-269.82 


39.35 


39.57 


87.26 


-9.73 


0.96 


73.00 


111.43 


-6.19 


18.52 


-321.23 


48.48 


61.29 


128.26 


-9.93 



TABLE IV. Cluster contributions to A18+UIX PNM E(p) in MeV. 



p 


T F 


(T)2B 


(T)mb 


AT 


<t> a >2B 


\V )2B 


(v s )mb 


\V ) MB 


5E2B 


0.02 


8.77 


3.89 


-1.66 


0.08 


-10.20 


0.11 


3.70 


-0.06 


-0.23 


0.04 


13.92 


5.15 


-1.77 


0.14 


-15.60 


0.38 


4.53 


-0.14 


-0.42 


0.08 


22.09 


8.43 


-2.25 


0.30 


-26.69 


1.27 


6.37 


-0.23 


-0.83 


0.12 


28.95 


11.97 


-2.70 


0.42 


-37.77 


2.51 


8.22 


-0.14 


-1.13 


0.16 


35.07 


15.90 


-4.00 


0.18 


-49.36 


3.87 


11.44 


0.04 


-1.26 


0.20 


40.69 


22.37 


0.42 


2.78 


-59.04 


10.65 


3.40 


4.87 


-1.29 


0.24 


45.95 


25.45 


1.15 


3.25 


-69.41 


13.03 


2.89 


6.32 


-1.65 


0.32 


55.67 


35.82 


2.49 


4.44 


-93.04 


18.04 


2.48 


11.04 


-2.45 


0.40 


64.60 


42.87 


4.34 


5.38 


-113.78 


25.07 


0.31 


16.94 


-3.20 


0.48 


72.95 


51.57 


6.35 


6.27 


-135.45 


32.61 


-2.18 


24.67 


-4.00 


0.56 


80.84 


61.73 


8.57 


7.07 


-157.94 


40.68 


-4.92 


34.25 


-4.75 


0.64 


88.37 


80.34 


12.00 


7.62 


-185.72 


53.35 


-8.56 


47.64 


-5.45 


0.80 


102.54 


101.95 


18.33 


7.90 


-230.35 


77.63 


-16.19 


74.42 


-6.50 


0.96 


115.80 


120.83 


21.91 


7.62 


-273.05 


99.09 


-22.17 


109.75 


-7.35 
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TABLE V. Contributions of pion-exchange and phenomenological parts of nuclear interactions 
to the E(p) in MeV. 





SNM 






SNM 






PNM 






PNM 








A18 






A18+ 


UIX 




A18 






A18+ 


UIX 




p 


K) 


<«§> 


(vfj) 


<«§> 








<«§> 


(vfj) 


<«§> 


<vgJE> 




0.04 


-11.4 


-6.4 


-11.8 


-6.4 


-0.4 


0.3 


-1.7 


-9.0 


-1.8 


-9.1 


0.1 


0.1 


0.08 


-19.2 


-12.1 


-19.1 


-12.1 


-0.8 


1.5 


-2.7 


-16.1 


-3.0 


-16.3 


0.3 


0.8 


0.12 






-27.5 


-16.0 


-2.1 


3.4 






-3.9 


-23.2 


0.6 


2.2 


0.16 


-32.6 


-20.9 


-35.3 


-19.8 


-3.6 


6.4 


-4.6 


-29.3 


-4.1 


-29.9 


1.2 


4.5 


0.20 






-42.0 


-22.8 


-5.5 


10.6 






-17.2 


-22.9 


-8.7 


10.1 


0.24 


-41.3 


-28.1 


-48.5 


-24.6 


-8.1 


15.9 


-5.8 


-41.1 


-19.1 


-28.0 


-10.1 


15.2 


n ^9 

U.oZ 


/I7 7 
-4 < . < 


-66. i 


-oy.o 


97 P. 


-16.6 


ou.y 


P. Q 
-0.6 


-DU.4 


-ZD. 1 


-OD.4 


17/1 
-1 1 .4 


oU.D 


0.40 


-54.7 


-36.7 


-72.6 


-17.4 


-38.4 


53.0 


-7.4 


-56.9 


-30.1 


-41.4 


-24.2 


50.9 


0.48 


-61.2 


-37.6 


-81.4 


-16.8 


-50.2 


80.3 


-9.6 


-63.4 


-35.7 


-44.7 


-34.1 


78.0 


0.56 


-65.7 


-37.0 


-91.5 


-12.0 


-65.2 


114.5 


-33.0 


-41.9 


-41.6 


-46.3 


-47.2 


112.7 


0.64 


-66.9 


-35.8 


-99.4 


-7.3 


-82.0 


155.9 


-39.2 


-42.2 


-52.2 


-41.1 


-76.9 


160.3 


0.80 


-67.1 


-40.5 


-112.1 


8.5 


-117.7 


260.2 


-42.3 


-34.4 


-63.3 


-31.2 


-116.3 


267.9 


0.96 


-70.9 


-19.5 


-121.2 


38.0 


-169.4 


397.7 


-46.2 


-11.9 


-70.9 


-15.5 


-155.0 


402.5 



TABLE VI. The E{p) of SNM in MeV. 



p 


A18 


A18+(5t; 


A18+UIX 


Al8+fe+UIX* 


corrected 


0.04 


-4.28 


-4.08 


-4.39 


-4.31 


-6.48 


0.08 


-8.72 


-8.07 


-8.06 


-7.97 


-12.13 


0.12 






-10.52 


-10.54 


-15.04 


0.16 


-14.59 


-12.54 


-11.85 


-12.16 


-16.00 


0.20 






-11.28 


-12.21 


-15.09 


0.24 


-17.61 


-13.69 


-8.99 


-10.89 


-12.88 


0.32 


-18.13 


-11.87 


0.84 


-4.21 


-5.03 


0.40 


-16.37 


-7.70 


12.23 


2.42 


2.13 


0.48 


-12.21 


-1.01 


32.18 


15.56 


15.46 


0.56 


-5.79 


8.16 


59.99 


34.42 


34.39 


0.64 


2.76 


19.54 


95.05 


58.36 


58.35 


0.80 


25.01 


45.24 


188.51 


121.25 


121.25 


0.96 


56.51 


82.63 


313.46 


204.02 


204.02 
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TABLE VII. The E(p) of PNM in MeV. 



p 


A18 


A18+Sv 


A18+UIX 


A18+<to+UIX* 


0.02 






4.35 


4.45 


0.04 


6.06 


6.32 


6.23 


6.45 


0.08 


8.53 


9.26 


9.21 


9.65 


0.12 






12.71 


13.29 


0.16 


12.33 


14.51 


17.38 


17.94 


0.20 






23.47 


22.92 


0.24 


16.69 


20.76 


28.85 


27.49 


0.32 


22.19 


28.59 


43.28 


38.82 


0.40 


29.41 


38.10 


63.79 


54.95 


0.48 


38.91 


50.35 


90.46 


75.13 


0.56 


49.08 


66.00 


123.93 


99.75 


0.64 


59.37 


81.15 


165.40 


127.58 


0.80 


88.27 


119.46 


273.37 


205.34 


0.96 


125.29 


167.02 


412.30 


305.87 



TABLE VIII. Contributions of relativistic boost interactions to the E(p) in MeV. 



p 


A18 

SNM 

{Sv) 2 B 


{5v)mb 


PNM 

(Sv) 2 B 


(Sv)mb 


A18+UIX 
SNM 

{Sv)2B 


{5v)mb 


PNM 

($v)2B 


(Sv)mb 


0.02 














0.08 


0.03 


0.04 


0.15 


0.05 


0.21 


0.05 


0.15 


0.05 


0.21 


0.05 


0.08 


0.45 


0.20 


0.59 


0.14 


0.44 


0.20 


0.58 


0.15 


0.12 










0.85 


0.40 


1.10 


0.29 


0.16 


1.32 


0.74 


1.71 


0.46 


1.33 


0.73 


1.72 


0.49 


0.20 










1.88 


1.09 


2.60 


0.58 


0.24 


2.42 


1.50 


3.16 


0.91 


2.49 


1.50 


3.42 


0.84 


0.32 


3.76 


2.49 


4.92 


1.48 


3.86 


2.51 


5.46 


1.41 


0.40 


5.18 


3.49 


6.64 


2.06 


6.04 


3.75 


7.73 


2.25 


0.48 


6.69 


4.52 


8.67 


2.77 


7.92 


5.17 


10.34 


3.20 


0.56 


8.35 


5.59 


12.28 


4.64 


9.96 


6.81 


13.27 


4.23 


0.64 


10.16 


6.62 


15.81 


5.97 


12.29 


8.69 


16.69 


4.81 


0.80 


12.61 


7.62 


22.12 


9.07 


17.12 


11.90 


23.83 


7.25 


0.96 


16.95 


9.17 


29.62 


12.11 


23.11 


14.61 


31.91 


10.59 
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TABLE IX. A18+UIX*: Contributions to (5v RE ) in MeV. 





SNM 








PNM 








p 


C2B 


W S {K) 


W S (F) 


w cch 


C2B 


W S {K) 


W S {F) 


w cch 


0.02 










0.04 


0.00 


0.02 


0.00 


0.04 


0.09 


0.00 


0.04 


0.00 


0.11 


0.00 


0.04 


0.00 


0.08 


0.24 


-0.02 


0.15 


0.00 


0.31 


-0.01 


0.12 


0.00 


0.12 


0.46 


-0.04 


0.31 


0.00 


0.58 


-0.03 


0.23 


0.00 


0.16 


0.73 


-0.06 


0.54 


0.00 


0.93 


-0.03 


0.39 


-0.01 


0.20 


1.03 


-0.09 


0.81 


0.00 


1.32 


-0.16 


0.50 


-0.03 


0.24 


1.36 


-0.14 


1.11 


0.00 


1.73 


-0.19 


0.68 


-0.03 


0.32 


2.11 


-0.22 


1.85 


0.00 


2.77 


-0.32 


1.14 


-0.06 


0.40 


3.25 


-0.70 


3.17 


-0.04 


3.90 


-0.40 


1.69 


-0.08 


0.48 


4.26 


-0.94 


4.35 


-0.04 


5.23 


-0.51 


2.34 


-0.09 


0.56 


5.36 


-1.15 


5.61 


-0.05 


6.71 


-0.66 


3.05 


-0.14 


0.64 


6.63 


-1.53 


7.20 


-0.05 


8.51 


-1.00 


3.68 


-0.24 


u.ou 


Q D9 


1 Qfi 


Q ^8 


-0.05 


12.14 


1 43 




n 39 


0.96 


12.44 


-3.58 


13.37 


-0.16 


16.23 


-1.87 


7.65 


-0.40 






TABLE X. 


A 1 O 1 T TT"V ^ 

A18+ULX : 


Contributions to (Sv 


) in MeV. 








SNM 








PNM 








P 


C2B 


W S {K) 


W S (F) 


w cch 


C2B 


W S {K) 


W S {F) 


w cch 


0.02 










0.03 


0.00 


0.01 


0.00 


0.04 


0.07 


0.00 


0.02 


0.00 


0.09 


0.00 


0.03 


0.00 


0.08 


0.20 


-0.02 


0.08 


0.00 


0.27 


-0.01 


0.06 


0.00 


0.12 


0.39 


-0.04 


0.17 


0.00 


0.51 


-0.01 


0.11 


-0.01 


0.16 


0.59 


-0.06 


0.32 


-0.01 


0.80 


-0.02 


0.18 


-0.02 


0.20 


0.85 


-0.07 


0.48 


-0.02 


1.28 


-0.12 


0.42 


-0.03 


0.24 


1.13 


-0.11 


0.65 


-0.03 


1.69 


-0.13 


0.55 


-0.04 


0.32 


1.76 


-0.16 


1.09 


-0.06 


2.68 


-0.20 


0.92 


-0.07 


0.40 


2.79 


-0.52 


1.96 


-0.14 


3.82 


-0.23 


1.37 


t\ -1-1 

-0.11 


0.48 


3.66 


-0.64 


2.67 


-0.23 


5.11 


-0.29 


1.93 


-0.16 


0.56 


4.60 


-0.72 


3.45 


-0.31 


6.55 


-0.37 


2.55 


-0.21 


0.64 


5.66 


-0.89 


4.39 


-0.43 


8.17 


-0.59 


3.19 


-0.25 


0.80 


7.79 


-0.96 


5.77 


-0.78 


11.69 


-0.85 


4.81 


-0.37 


0.96 


10.67 


-1.49 


7.80 


-1.34 


15.68 


-1.06 


6.82 


-0.55 
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TABLE XI. Maximum gravitational masses, in M Q , and moments of inertia /, in M & km 2 , for 
stars with beta-stable matter and with PNM; with incompressible (INC) matter at p > pj (fm -3 ) 
for the Al8+(5u+UIX* model; and with mixed nuclear-quark matter (NM+QM) phase (bag con- 
stant B in MeV fin" 3 ). 



NM 


Max. mass 


Max. mass 


Max. I 


pi of INC 


Max. mass 


Max. I 


models 


beta-stable 


PNM 


beta-stable 


Models 


beta-stable 


beta-stable 


Al8+<5u+UIX* 


2.20 


2.21 


115. 


0.32 


2.92 


261. 


A18+UIX 


2.38 


2.39 


143. 


0.48 


2.46 


157. 


A18+(5t; 


1.80 


1.81 


67. 


0.64 


2.26 


123. 


A18 


1.67 


1.68 


55. 


0.86 


2.19 


115. 


FPS 


1.80 




73. 








NM+QM 


B 


Max. mass 


Max. I 


B 


Max. mass 


Max. I 


models 




beta-stable 


beta-stable 




beta-stable 


beta-stable 


Al8+<5u+UIX* 


122 


1.91 


96. 


200 


2.02 


107. 


A18+&; 


122 


1.74 


66. 


200 


1.76 


67. 



TABLE XII. Parameter values for effective Hamiltonian for different models. The parameters 
P3 = 89.8 MeV fm 5 , p 4 = 0.457 fm 3 and p$ = —59.0 MeV fm 5 (see Eq. ([Ad])) are common to 
all four models. The dimensions of the parameters, involving MeV and/or powers of fm, can be 
worked out from Eqs. flA~l]) ; ( [All) or (|A4|). 



model 


Pi 


P2 


P% 


Pi 


P8 


P9 


Pio 


PU 


P12 


P13 


Al8+<fe+UIX* 


337.2 


-382. 


-19.1 


214.6 


-384. 


6.4 


69. 


-33. 


0.35 


0. 


A18+UIX 


328.8 


-404.6 


-34. 


217.5 


-385.6 


6.35 


25.4 


0. 


0.47 


-0.9 


A18+5v 


281.0 


-151.1 


-10.6 


210.1 


-158. 


5.88 


58.8 


-15. 


-0.2 


-0.9 


A18 


297.6 


-134.6 


-15.9 


215.0 


-116.5 


6.42 


51. 


-35. 


-0.2 


0. 


model 


PU 


P15 


P16 


Pn 


P18 


P19 


P20 


P21 






Al8+<fe+UIX* 


0. 


287. 


-1.54 


175.0 


-1.45 


0.32 


0.195 


0. 






A18+UIX 


-452. 


217.1 


-1.0 


100.3 


-1.19 


0.32 


0.2 


-275. 
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P [fm" 3 ] 
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